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1  INTRODUCTION 


1.1  GENERAL 

1.1.1  Contract/ Scope 

This  work  was  accomplished  under  Contract  Numnc-t  DAADOT-'iu-C- 

with  the  U.S.  Army  White  Sands  Missile  Range.  The  contract  became 

effective  on  7  February  1980  and  constituted  less  than  one  half  man  yea 

of  -echnica!  effort.  The  Contracting  Officer's  Technical  Representative  i 

Mr.  Ronald  Cionco 
Atmospheric  Sciences  Laoora’ory 
Attn:  DELAS-BE-C 

White  Sands  Missile  Range,  NM  8f  002 


7 


1.1.2  Overview 


Environmental  condition?  often  dictate  effectivener  s  :  Army  -  • 

and  defensive  weapon  systems.  Such  things  as  terrain ,  vccetaticr. ,  i-:  • 

weather  determine  in  what  direction,  how  far,  and  what  decree  of  clar:'-  : 
target  can  be  acquired,  isolated,  and  distinguished  from  its  surroundin'.:  . 
For  a  particular  season  and  theater  of  operation,  weather  offers  the  cr¬ 
est  day  to  day  change  that  influences  weapon  system  performance  . 
are  being  conducted  to  develop  techniques  to  determine  "seeabihty"  : 
obscuration  from  natural  as  well  as  battiefield  induced  eontamina’es  .  J 
siderable  attention  is  being  directed  towards  electro-optical  (EO)  energy 
absorption  due  to  gaseous  molecules,  scattering  due  to  haze,  s:no-;e,  f:o, 
cloud,  rain,  snow  or  hail,  and  beam  wandering  due  to  small  scale  tur: 
lence  and  temperature  fluctuations.  Obscuration  due  to  scattering,  howo\  o 
depends  greatly  upon  the  energy  wavelength  of  *he  viewing  sensor  relati.-e 
to  type,  number  and  size  distribution  of  atmospheric  litno-  and  hydro- 
meteors  , 

The  scattering  process  itself  behaves  nonlinearly  and,  by  necessity, 
obscuration  models  must  consider  such  meteorological  microphysical  fea¬ 
tures  as  number  concentration  and  size  distribution  of  particles.  It  is 
unreasonable  to  expect  observations  of  these  microphysical  features  within 
a  battlefield  environment.  Thus  it  becomes  necessary  to  work  within  the 
scope  of  routinely  available  weather  information  to  infer  the  microfea- 
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tures  and  to  define  the  actual  scattering  medium .  Geo-Atmospherics  Cornu- 
ration  has  developed  weather  models  to  depict  hour  oy  hour  variations  o: 
scattering  media  important  to  Army  weapon  system  performance. 

Better  instruments  and  experimental  data  ga’hermg  programs  are  nec¬ 
essary  to  establish  reliable  data  bases  of  atmospheric  microohysical  fea’urc; 
of  liquid  and  solid  particles  associated  with  the  wide  range  of  possible 
global  weather.  These  data  must  have  the  accuracy  demanded  for  basic 
studies  on  their  impact  on  design  and  operation  of  EO  systems.  To  ga’e, 
a  number  of  static  models  have  been  developed  t u  reprr  son*  atmospben  - 
particles  for  a  particular  type  of  environment.  for  example,  aerosol  mo  :•  1. 
have  been  generated  to  represent  "typical"  conditions  expected  to  be  f-_r.  i 
m  Continental  or  maritime  regions  or  at  different  levels  m  T.e  atmos  p.-mr1'  . 
Some  of  these  models  are: 

A.  Rural  Model  - 

B.  Urban  Model  - 

C.  Maritime  Model  - 


I 

which  is  to  represent  the  natural 
midlatitude  environment  found  in  the 
country  or  in  clean  air  urban  region; 
found  after  a  cold  front  passage. 

which  is  to  represent  the  industrial  [ 

aerosols  found  in  cities  and  in  rural 
regions  experiencing  air  stagnation 
and  subsequent  air  pollution  build-uo. 

which  is  to  represent  the  open  ocean 

regions  at  least  300  km  offshore  wrh  j 

a  moderate  surface  3ir  wind  speed. 
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D.  Tropospheric  Model  -  which  is  to  normally  represent  T.e 

turbulent  layer  of  the  atmosphere  be¬ 
tween  the  top  of  the  boundary  layer 
(about  2  km)  and  the  tropopause  (u 
18  km,  depending  upon  season  and,  t 
a  greater  extent,  latitude)  as  well  as 
those  special  boundary  layer  case? 
when  the  surface  air  layer  over  land 
is  calm  and  clean  (meteorological 
range  greater  than  40  km). 

E.  High  Level  Models  -  which  is  to  represent  aerosol  distri¬ 

butions  within  the  stable  stratosphere 
(from  the  tropopause  to  about  2  0  km' 
and  the  mesosphere  (from  the  top 
the  stratosphere  to  about  80  km). 

Similar  types  of  static  models  have  been  developed  for  fog,  cloud, 
and  rain  conditions.  The  atmosphere  is,  however,  a  fluid  in  motion  that 
produces  continuous  rather  than  discrete  spectrums. 

What  we  have  done  here  is  to  take  the  first  step  in  developing  dy¬ 
namic  models  that  relate  particle  type,  concentration,  size  distribution,  and 
their  vertical  variation  to  observed  dynamics  and  thermodynamics  of  the 
atmosphere  that  exist  at  any  instant  in  time.  A  ground  rule  established  was 
that  only  routinely  available  standard  surface  and  upper  air  observations 
and/or  forecasts  would  be  available  for  use  in  a  real-world  battlefield  en¬ 
vironment.  As  such,  a  number  of  simplifying  assumptions  had  to  be  made, 
to  use  meso  and  synoptic  scale  data  to  depict  microscale  features.  Not 
all  gaps  were  closed  and  much  remains  to  be  done  but  significant  advance? 


were  made  in  dynamically  modeling  litho-  and  hydro-meters  important 
operational  performance  of  Electro  -  Optical  Systems. 


11 


1.2  MODELING  EQUATIONS 
1.2.1  General 

There  is  no  question  that  adverse  weather  is  detrimental  to  electro¬ 
magnetic  energy  propagation.  Theoretical  and  application  techniques  have- 
been  developed  to  simulate  atmospheric  radiative  transfer  effects  to  derive 
obscuration  expected  for  any  sensor.  What  is  required  is  a  qualitative 
description  of  the  scattering  medium  that  is  needed  in  order  to  be  able 
to  derive  obscuration.  The  problem  is  how  to  map  three  dimensional  dis¬ 
tributions  of  microphysical  scattering  media  given  only  routinely  available 
surface  and  upper  air  meteorological  information. 

Before  dicussing  the  techniques  used  to  develop  the  dynamical  models 
for  fog,  cloud,  and  rain  conditions,  it  is  worthwhile  to  consider  rhe  dif¬ 
ferent  characteristics  associated  with  the  two  most  popular  analytical 
methods  used  to  depict  particle  size  distributions,  namely  the  power  law 
and  the  exponential  (modified  gamma)  equations.  The  equations  will  be 
presented  and  discussed  to  show  shapes  and  characteristics  responsive 
to  the  wide  variety  of  atmospheric  particles.  For  the  purposes  of  this 
study  the  particle  sizes  vary  by  six  orders  of  magnitude.  Each  type  of 
particle  has  a  typical  mode  radius  and  range  of  values,  as  shown  in 
Table  I.  Although  size  distributions  vary  within  each  category  which  must 
be  considered,  it  is  equally  important  to  develop  modeling  schemes  that 
account  for  the  extreme  variations  among  different  meteorological  events. 
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TABLE  I  Typical  Atmospheric  Particle  Sizes 


i 

! 

► 

I 

i 

i 

! 


i 


Particle  Type 

Mode  Radius 

micrometers 

Radius  Rang 

micrometers 

Continental  Aerosols 

0.  02 

.002  -  20 

Maritime  Aerosols 

0.2 

.02  -  20 

New  Fog  (  1  hr  ) 

4 

i  -  40 

Old  (Evolved)  Fog  (  2  hr  ) 

10 

1  -  °0 

Fair  Weather  Cumulus 

3.5 

1  -  1  j 

Stratus/Stratocumulus 

4 

1  -  15 

Altos  tratus 

5 

1  -  1  3 

Nimbostratus 

5 

1  -40 

Cumulonimbus 

6 

1  -  100 

Mist  (0.05  mm/hr) 

75 

6  -  5  5  0 

Drizzle  (0.25  mm/hr) 

150 

6  -  7  50 

Light  Rain  (1  mm/hr) 

175 

6  -  1250 

Moderate  Rain  (4  mm/hr) 

200 

o  -  1750 

Heavy  Rain  (16  mm/hr) 

225 

6  -  225'"' 
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1.2.2  Power  Law  Si2e  Distribution 


One  of  the  most  popular  methods  in  use  today  to  represent  solid 
aerosols  in  the  atmosphere  was  developed  by  Junqe  in  1953.  Prom  many 
of  his  aerosol  measurement  studies  he  found  the  concentration  pea'<ed  at 
a  radius  of  about  0.02  n m  hut  remained  linearly  distributed  between  about 
0.1  to  10  fim,  when  plotted  on  a  log  log  scale.  This  linear  range  covers 
the  most  optically  interesting  phenomena  in  a  hazy  atmosphere.  The 
smaller  aerosols  in  this  range  provide  a  bluish  tinge  to  thin  haze,  those 
around  0.3  /urn  are  most  important  in  determining  what  we  call  visibility 
in  haze,  and  the  larger  aerosols,  which  are  fewer  in  number,  can  provide 
beautiful  red  sunsets.  Therefore,  Junge's  power  law  representation  of 
aerosol  size  distribution  and  concentration  is  very  appropriate  in  defining 
E0  responses  to  haze  and  it  is  simple  to  derive  from  measured  data  an 
simple  to  use. 

Junge's  power  -  law  size  distribution  function  is: 

"W  -  c  r  "v 

where  n(r)  is  the  number  of  particles  per  unit  interval  of  radius  and  per 
unit  volume,  C  is  the  surface  air  particle  concentration  factor,  r  is  'he 
aerosol  radius  in  fim,  and  the  exponent  v  defines  the  slope  of  the  distri¬ 
bution  curve. 
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Often  the  nonlogarithmic  form  of  equation  (1)  is  desired.  Since 

d  log  r  =  0.434  d  In  r  and  d  In  r  =  ~ 

then  equation  (1)  can  be  rewritten  as: 

n  (r)  =-^-  =  0.434  C  r  _(  v  *  1  ]  (2) 

We  will  use  this  form  of  the  power  -  law  equation.  Therefore  the  particl 
concentration  (  N  )  say  for  all  particles  equal  to  or  greater  than  radius  r, 
becomes  the  integral  of  equation  (2  )  as 

CL 

N  =  0.434  C  f  r  '  {  v  +  1  )  d  r  (3) 

rl 

The  units  for  droplet  or  aerosol  concentrations,  n  (r) ,  are  cm-^  //m  . 

Therefore  the  units  for  total  particle  concentrations  above  a  given  radius 
are  cm"^, 
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Fig.  1.  Power  Law  Size  Distributions 


6 


Figure  1  shows  how  the  power  law  size  distribution  given  in  Eq.  (2) 
varies  as  a  function  of  the  exponent  v  and  the  factor  C.  Plotted  in  Fig.  1 
is  the  aerosol  concentration  as  a  function  of  haze  particle  size  radius. 
These  plots  are  made  for  continental  air  near  the  earth’s  surface  assuming 
two  different  slope  factors  v  and  two  visibilities,  which  effect  the  con¬ 
centration  factor  C.  Later  we  will  describe  the  rational  and  equations  wo 
developed  to  specify  C  and  v  according  to  changing  meteorological  features. 
For  now,  we  want  to  demonstrate  how  the  power  law  distribution  changes 
with  a  fixed  C  and  variable  v  and  for  a  fixed  v  with  variable  C.  Fig.  1 
shows  two  visibilities,  1  and  10  km,  which  translates  to  concentration 
factors  C  of  30  and  3,  respectively.  The  exponents  shown  correspond  to 
v  values  of  3  and  2. 

Notice  in  Fig.  1  for  a  constant  concentration  factor,  visibility  1  km, 
as  the  exponent  slope  factor  v  decreases  in  value  the  aerosol  size  distri¬ 
bution  curve  becomes  more  shallow  resulting  in  fewer  numbers  of  smaller 
particles  and  larger  concentrations  of  larger  size  particles.  Typical  values 
of  v  are  2  and  3  for  fog  and  haze  particles,  respectively,  and  are  shown 
plotted  In  Fig.  1  for  1  km  visibility  conditions.  It  is  not  unusual  to 
have  both  haze  and  fog  present  during  1  km  visibility  conditions.  Fig.  1 
also  shows  that  curves  for  a  fixed  value  of  v  are  parallel  but  displaced 
vertically  according  to  the  value  of  the  concentration  factor  C  which 
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represents  the  number  of  particles  restricting  the  visibility.  Thus  by  r 
lating  the  value  of  v  and  C  to  meteorological  events  it  is  possible  to 
derive  corresponding  particle  size  distributions. 


18 


p 


1.2.3  Exponential  or  Gamma  Distribution 

One  disadvantage  of  the  above  described  power  law  distribution  func¬ 
tion  is  that  the  curve  is  linear  on  a  log  log  plot.  This  means  that  if  a 
preferred  concentration  of  particle  sizes  occurs  then  a  peak  would  exist  in 
the  distribution  and  may  not  be  well  represented  by  the  power  law  approacr. 
Haze,  fog,  cloud,  and  rain  particles  often  have  a  Pearson  type  111  or  log¬ 
normal  distribution  where  the  number  of  particles  increases  rapidly  as  parti¬ 
cle  radius  increases  and  reaches  a  peak  value  after  which  there  is  a 
slower  decay  or  trailing  off  of  particles  at  increasingly  larger  sizes.  The 
exponential  or  gamma  distribution  can  handle  these  types  of  distributions 
but  at  the  expense  of  greater  complexity  than  the  power  law  distribution. 

A  popular  form  of  the  exponential  function  is  given  by  Deirmendjiar. 
as 


n  (r)  =  d  N  /  d  r  =  A  r°^xp  (  -  B  r  y  )  (4) 

where  n  (r)  units  are  cm  ^  /ira'^,  with  particle  radius  r  in  micrometers, 
ans  A,  B,a,  and  7  being  positive  constants.  Taking  the  derivative  of 
Eq.  (4)  and  settling  it  equal  to  zero  gives  the  mode  radius,  r  ,  that  rep¬ 
resents  the  radius  of  maximum  concentration  occuring  at  the  peak  of  the 
size  distribution  curve  and  is  expressed  analytically  as 
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r£  =  <*/  Br  (5) 

The  concentration,  N,  or  total  number  of  particles  per  uni*  volume 
can  be  derived  by  taking  the  integral  of  Eq.  (4)  between  the  limits  of 
zero  and  infinity  to  obtain 

N  =  Ay-1  b  1  ]  /y  r(-SLLl)  (6) 

where  /"represents  the  gamma  function.  From  this  we  can  see  that  the 
coefficient  A  is  proportional  to  the  visibility  or  number  concentration  of 
atmospheric  particles  and  is  given  by 

a  =  n  y  b  (a't’  1/y)  /  r(a+  1  / y )  (7) 

Most  often  y  is  set  equal  to  1  for  computation  simplicity  and  to  reduce 
Eq.  (6)  so  that  the  gamma  function  is  finite  if  a  -*■  1  >  0  and  if  ct  - 
is  an  integer  then 


p{ct  +  1  )  =  a  !  (8) 

Thus  with  the  above  simplifications  of  Y  =  1  and  a  +  1  an  integer,  Equations 
(4),  (5),  (6),  and  (7)  can  be  rewritten 

n  (r)  =  d  N  /  d  r  5  A  r"  exp  (  -  B  r  )  (9) 


rc  =  a  /  B 

(10) 

N  =  A  a  !  B  1  } 

(ID 

J 


i 
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and 


A  =  (  N  /  a  !  )  B 


(  a  +  1  ) 


(12) 


Normally  what  is  done  is  to  select  constant  values  for  A,  F,  or,  N , 
and  rc  to  represent  static  conditions  for  a  given  fog,  cloud,  or  tain  type- 
in  order  to  specify  particular  models  to  represent  "average"  particle  size 
conditions.  We  are  interested  here  in  developing  dynamical  models  of 
particle  size  distribution  by  turning  the  above  so  called  "constants"  into 
variables  that  are  related  to  the  type  and  changes  of  actual  weather  ob¬ 
servations,  Therefore  let  us  now  attempt  to  translate  these  "constants" 
in  terns  of  meteorological  variables. 

The  linear  mean  radius,  r  ,  is  the  sum  of  all  ’he  drorle-  radii  di¬ 
vided  by  the  total  number  of  droplets  so 


rm 

f  r  n  (r)  d  r 

N  J0 

(12) 

= 

1  A  ( a  +  1  )  ! 

N  B  (a  +  2  ) 

(14) 

so 

rm 

(a-  1  )  / B 

(1  5) 

therefore 

rm 

(on-  l  )  rc  /a 

(lb) 

By  definition  the 

fog  or  cloud  liquid  water  content,  W^, 

gives  1 1 

total 

mass  concentration  of  liquid  water  by 

(  g  m-3  ) 

=  10  "6  (  4"/  3  )  P  /*r3  n(r)  d  r 

Jo 

(17) 
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where  P  is  the  density  of  water  in  g  cm  and  r  is  m  microme'ers. 
Expanding  Eq.  (17)  gives 

WL  =  10"6  (  4  n/  3  )  Pw  A  (a-  3  )  !  /  B  4  ]  (1  -i) 

Meteorologically  speaking  tlie  density  of  water  is  essentially  one  c  cr.*', 
the  liquid  water  content  can  be  inferred  from  visibility  observations  lr.  fog, 
the  particle  displacement  along  a  moist  adiabat  in  clouds,  and  the  in¬ 
tensity  for  rainfall,  and  the  coefficient  B  is  defined  by  the  mode  radius  as 
shown  in  Eq.  (10).  Thus  we  can  solve  for  the  coefficient  A  in  Eq .  (IB) 
and  obtain 

A  =  10+6  (4n/3)~l  (a/  rc  )(a+  4  *  WL  /  (a-  3  )  !  (19) 

Qoud  and  fog  dropsize  observations  show  that  the  mode  radius  rc  varies 
with  such  things  as  cloud  or  fog  type,  age,  and  height  above  cloud  or  fog 
Dase.  Assuming  we  can  reasonably  approximate  the  mode  radius  and  liquid 
water  content  then  we  are  left  with  only  the  shape  factor  a  as  a  variaele 

needed  before  we  can  specify  particle  size  distributions  given  by  Eq.  (9). 

We  have  made  a  number  of  calculations  to  show  how  the  particle 
size  distribution  changes  according  to  individual  changes  in  liquid  water 
content,  mode  radius,  and  curve  shape  factor.  First,  Fig.  2  shows  two 
different  ways  used  to  plot  particle  size  distributions .  The  top  log- log 
plot  is  a  computer  output  solution  of  the  above  gamma  equation  for  five 
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RADIUS  (micromct«r»l 


Two  Ways  to  Plot  Cloud  Particle  Size  Distributions 


different  cloud  types.  The  bottom  linear  plot  of  observed  particles  for 
altostratus  (As),  stratocumulus  (Sc),  and  nimbostratus  (Ns)  clouds  crows 
the  familiar  Pearson  type  III  distribution.  There  is  no  preferrcj  way  to 
present  such  data.  We  will  use  mostly  the  log-loq  method  because  sma 
ler  concentrations  of  particles  can  be  plotted  with  greater  ease. 
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Figures  3,  4,  and  5  show  how  the  particle  size  distribution  giver, 
by  the  exponential  function  varies  when  only  one  variable  is  allowed  to 
change.  Plots  in  Fig.  3  depict  two  orders  of  magnitude  change  in  liquid 
water  content  while  the  mode  radius,  rc  =  4/im,  and  shape  factor,  or  =  4, 
remain  constant.  What  is  most  striking  in  Fig.  3  is  that  the  maximum, 
density  d  N  /  d  r  varies  from  about  4.5  to  450  droplets  per  uni:  volume  and 
radius  interval  (  cm'^,  )  which  represents  significant  differences  m 

density  of  or  visibility  in  clouds  or  fogs.  Thus  if  we  can  infer  the  mode 
radius  and  shape  factor  for  a  specific  type  of  fog  or  cloud  then  we  can 
infer  the  number  density  of  the  droplets  by  obtaining  a  measure  of  the 
liquid  water  content. 

We  made  calculations  to  show  how  variations  in  the  shape  factor 
coefficient  a  could  be  used  to  represent  narrow  or  broad  distributions  m 
particle  sizes.  Fig.  4  shows,  for  a  fixed  liquid  water  content  and  mode 
radius,  that  as  one  selects  higher  and  higher  values  for  a  one  obtains 
sharper  and  sharper  curves  that  represent  a  more  narrow  distribution  of 
particle  sizes.  Also  note  that  as  the  drop  size  distributions  become  broad¬ 
er,  the  number  density  decreases  for  the  most  frequently  occurring  particles 
while  increasing  for  the  small  and  large  particles.  Therefore,  by  knowing 
something  about  the  variance  or  distribution  of  particle  sizes  associated 
with  particular  fogs  or  clouds  it  should  be  possible  to  select  an  appropriate 
shape  factor. 
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RADIUS  ( micrometers ) 


Exponential  Function  Size  Distributions  For  Three  Liquid  Water 
Contents,  WL=  0.02,  0.2,  and  2.0  g  m"3,  With  Constant  Mode 
Radius,  rc*  4  micrometers,  and  Shape  Factor,  a=4. 


RADIUS  (micrometers) 


Exponential  Function,  Size  Distributions  Tor  Three  Shape  factors, 
0 1-  2,  4,  and  (> ,  With  Constant  Liquid  Water  Content,  WL~0.2gm 
and  Mode  Radius,  r„  =  4  micrometers. 


Exponential 

rc  ”  ^  ® ' 

Content,  W. 


RADIUS  (micrometers) 


Function  Size  Distributions  For  Throe  Mode  Radii 
and  10  micrometers,  With  Constant  Liquid  Water 


=  0,2  g  m 


-3 


and  Shape  factor,  a  -  4. 


Lastly,  we  calculated  droplet  density  for  three  environments  wr.erc- 
the  most  frequently  occurring  particle  si2e  was  4,  6,  and  10  micrometers. 
During  these  computations,  the  liquid  water  content  and  shape  factor  were 
held  constant.  The  results,  Fig.  5,  show  the  peak  curve  displacement 
corresponds  with  the  selected  mode  radii  and  the  peak  number  density  of 
particles  decreases  as  mode  radius  increases.  Often  a  particular  mode 
radius  is  representative  of  a  particular  type  of  fog  or  cloud.  Thus  by  hav¬ 
ing  information  or  observations  on  the  liquid  water  content,  the  radius  of 
the  most  frequently  found  droplet,  and  the  variance  or  breath  of  particle 
sizes,  we  can  select  variables  in  the  exponential  or  gamma  function  equa¬ 
tion  to  depict  a  wide  range  of  unimodal  particle  size  distributions. 

In  the  following  report  we  will  be  using  both  the  power  law  and 
exponential  size  distribution  concepts  described  above  to  develop  lit'no- 
and  hydro-meteor  models  that  can  be  specified  given  only  routinely  avail¬ 
able  weather  observations. 


2  DYNAMIC  MODELS  OF  PARTICLE  SIZES 


2.1  LITHOMETEORS 

2.1.1  Continential 

2. 1.1.1  Ground  Level  Variations 

The  underlying  surface  (land,  ocean,  vegetative,  barren),  soil  type 
(humus,  clay,  rocky),  soil  condition  (wet,  dry,  compactei,  loose),  atmo¬ 
spheric  stability  (convective,  stable),  and  wind  speed  (soil  erosion,  vertical 
and  horizontal  advection)  are  important  in  determining  type,  concentration, 
and  size  distribution  of  lithometer  particles  as  a  function  of  space,  time, 
and  height  in  the  atmosphere.  To  further  complicate  matters  some  dust 
grains  are  nonhydroscopic  while  others,  especially  salt  containing  particle? 
are  highly  hydroscopic.  These  hydroscopic  particles  include 
many  well  known  things  as  salt  from  the  sea,  tars  and  resins  from  organic 
plants,  smoke  from  industrial  and  battlefield  fires,  and  the  usual  sulfate 
and  nitrite  chemicals  from  organic,  combustion,  and  photochemica  1  processe 
The  thing  is,  hydroscopic  particles  are  not  anly  important  as  light  scat¬ 
tered  but  also  provide  particle  size  changes  as  the  ambient  relative  humidity 
approaches  saturation  and  also  provide  the  most  active  condensation  nuclei 
needed  in  the  cloud  physics  process  of  cloud/fog  droplet  and  precipitation 
formation.  In  fact  studies  have  been  conducted  to  show  light  attenuation 
or  visual  range  restriction  as  relative  humidity  increases,  Fig.  b.  The 
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important  fact  to  us  is  that  visual  range  or  visibility  is  nearly  linearly 


related  to  variations  in  relative  humidity  greater  than  about  seventy  per 
cent.  That  is  to  say  that  we 


flQ  •  Variation  In  Visual  Range  With  Relative  Humidity. 


will  consider  variations  in  the  particle  size  as  a  function  of  relative  hu¬ 
midity  to  be  incorporated  in  meteorological  observations  of  visual  range 
or  visibility. 

Dust  particles,  lithometeors ,  or  haze  are  usually  contained  within 
the  radius  size  range  of  0.01  to  10  micrometers  with  a  peak  or  most  fre¬ 
quent  particle  radius  of  about  0.02  micrometers.  There  are  two  points  of 
most  interest  to  this  study.  First,  most  Army  electro-optical  systems  are 
degraded  or  adversely  effected  by  particle  sizes  equal  to  or  greater  'han 
0.1  micrometer  in  radius.  Second,  particle  size  observations  at  the  earth’s 
surface,  at  different  geographical  locations  (Baltimore,  Minneapolis,  and 
Seattle)  representing  east  -  west  coast  and  mid-sections  of  U.S.A.,  and 
day  and  night  aerosol  observations  (Los  Angeles)  all  show  distributions 


having  a  r-4  dependence  as  given  in  Eq.  (2).  That  is  to  say,  the  con¬ 
tinental  aerosol  size  distribution  at  the  earth’s  surface  follows  a  power 
law  curve  whose  exponent  is  -  4.  The  exponent  in  Eq.  (2)  is  -(  v  -  1), 
therefore,  a  value  of  v  =  3  is  typical  of  ground  level  haze  aerosols.  The 
number  density  or  vertical  displacement  of  this  curve,  however,  can  vary 
considerably  depending  upon  the  number  of  particles  suspended  m  'he  air 
at  a  given  instant  in  time,  which  in  turn  is  related  to  meteorological  vis¬ 
ibility. 

Before  Eq.  (2)  can  be  solved  to  derive  the  number  density  of  haze 
aerosols  it  is  necessary  to  obtain  a  measure  of  the  particle  concentration 
factor  C.  McCartney  has  shown  that  this  concentration  factor  is  directly 
proportional  to  the  backscatter  coefficient  which  in  turn  is  propoitional  to 
the  meteorological  range,  R^ ,  at  optical  wavelengths.  We  have  plotted 
his  values  of  the  surface  air  particle  concentra tion  factor  (C)  as  a  function 
of  meteorological  range(  R  )  in  kilometers,  see  Fig.  7.  McCartney  pre¬ 
sented  data  points  for  v  in  Eq.  (2)  equal  to  2.5  and  4.  From  Fig.  7  u 
can  be  seen  that  our  derived  equation 

C  =  30/  Rm  (20) 

represents  a  reasonable  fit  to  these  data.  More  importantly,  Eq.(20)  al¬ 
lows  us  to  use  visibility  observations  reported  over  standard  meteorological 
networks  to  infer  variations  in  the  number  density  of  haze  par’icles 
as  weather  conditions  vary. 


2. 1.1.2  Upper  Level  Variations 

Aerosol  variations  with  respect  to  geographical  location  and  meteoro¬ 
logical  conditions  decrease  with  height  and  often  are  considered  to  be 
non-existent  near  the  troposphere,  about  10  km.  A  bimodal  distribution  of 
aerosols  is  sometimes  observed  to  increase  with  height  above  the  earth's 
surface.  In  a  static  time  dependent  atmosphere  it  is  not  unusual  for 
small  particles  to  coagulate  and  form  an  ever-increasing  size  of  lower 
concentration  particles.  There  are,  however,  atmospheric  regions  where 

other  forces  retard  this  action  and  the  maintenence  or  produc'ion  o:  small 

* 

•.articles  remains  high.  Such  appears  to  be  the  case  in  or  around  the 
tropopause  height.  This  effect  is  shown  in  part  by  observations  taken  bv 
B 1  if  ford  (1970)  over  Death  Valley,  as  shown  in  Fig.  9,  for  three  different 
heights,  namely  3.5,  6,  and  9km.  This  figure  also  shows  the  gravita'iona 
settling  effect  on  large  particles  so  that  a  fewer  number  of  larger  parti¬ 
cles  exist  at  higher  altitudes. 

A  great  deal  of  dispersion  exists  among  investigators  in  terms  of 
what  the  vertical  distribution  is  for  atmospheric  aerosol  particle  concen¬ 
trations.  This  is  clearly  illustrated  by  Fig.  9.  For  example,  the  spread 
at  12km  shows  a  dispersion  of  more  than  an  order  of  magnitude.  A  typical 
"clear"  atmospheric  profile  of  aerosol  concentration  as  a  function  of  heigh' 
is  presented  in  Fig.  9  as  plotted  values  of  "C" ,  which  correspond  to  a 
clear  atmosphere  having  a  23km  range  of  visibility.  Superimposed  upon 
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this  "clear"  plot  are  a  series  of  "H  "  values  representing  a  hazy  atmo¬ 
sphere  whose  influence  extends  from  the  ground  surface  to  5  km,  above 

which  the  "hazy"  and  clear  atmosphere  have  the  same  aerosol  number  con¬ 
centrations.  For  this  case,  the  hazy  atmosphere  correspond-  with  a  grourv: 

•  I 

level  horizontal  visibility  of  5km.  Also  plotted  for  comparison  is  Toon's 
and  Pollock's  data  (  T  ,  5  )  for  aerosol  number  concentrations  for  particles 
equal  to  or  greater  than  0.15  micrometers  in  radius  that  are  found  be¬ 
tween  the  surface  and  12km.  Elterman's  early  data  in  1964  resulted  m  a 
straight  line  in  Fig.  9  extending  from  the  surface  to  10km,  cut  more  recent 
•cat a  in  1968  with  better  instrumentation  displayed  a  break  at  4  km.  with  a 
small  change  or  nearly  constant  value  of  particle  concentration  wrf.  heigr/, 
see  Fig.  9.  As  shown  in  Fig.  9,  the  hazy  and  clear  models  as  veil  as 
the  Elterman  and  Toon's  data  plots  include  those  particle  sizes  that  arc 
most  important  in  either  remotely  probing  the  atmosphere  or  studying  op¬ 
tical  or  infrared  radiational  characteristics  of  the  atmosphere.  What  is  no- 
shown,  however,  is  how  these  particle  size  distributions  change  with  heigh 
and  changing  weather  environments. 

We  were  looking  for  a  measure  of  the  change  in  particle  size  dis¬ 
tribution  as  altitude  increases  above  ground  level.  Bhfford’s  balloon  borne 
observations  over  Death  Valley,  Fig.  8,  showed  that  the  particle  density 
not  only  decreased  with  height  but  so  did  the  slope  of  the  particle  dis- 
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tribution  curve.  More  recent  aircraft  data  were  collectej  over  V/,- 
Europe  by  Cress,  see  fig.  10,  which  showed  similar  result  •;  r  ;• 

general  the  observed  distributions  near  the  surface  had  a  r”'1  while 
at  about  5km  had  a  r-^  slope.  Cress  plotted  the  exponent  of  r  as  w 
Junge  Slope  as  a  function  of  height  and  for  spring.  Simmer,  ar.ti  fail. 
obtained  a  linear  fit  to  the  fall  data.  Fig.  11,  and  derived 

v  =  3  -  /.  /  2  (2! ) 

where  the  Junge  Slope  =  -  (v  *  1  )  as  given  in  the  power  -  law  Eg.  (2), 
and  Z  is  the  height  about  ground  surface  in  km..  The  airborne  equipment 
sampled  aerosol  sizes  from  0.2  to  b  micrometers.  Observations  were 
taken  under  a  wide  range  of  meteorological  conditions,  including  clear  anti 
overcast  skies  and  visibilities  varying  from  excellent  down  *o  3km. 
scatter  of  data  points  exists  at  all  heights  uud  should  be  expected  rue  • 
physical  variations  previously  discussed,  e.q.  underlying  surface,  atmo¬ 
spheric  stability,  wind,  etc.  Cress  summer  Junge  Slope  observations  have 
about  the  same  slope  found  in  the  fall  data  but  the  surface  value  for  the 
Junge  Slope  is  -5.  This  implies  that  the  warmer  and  dryer  underlying 
surface  and  more  unstable  boundary  layer  during  summertime  generates  a 
larger  number  of  smaller  aerosols.  The  spring  observations,  Fig.  12,  had 
a  curvalinear  distribution  of  Junge  Slope  with  height  and  no  one  linear 
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Aircraft  Aerosol  Observations  at  Different  Heights 


JUNGE  SLOPE 

(for  particles  larger  than  4  micrometer  radius) 


Fig. 12 


Spring  Variation  of  Junge  Slope  With  Height 


curve  covers  the  entire  range  of  data.  For  the  fail  and  spring  data  i 


linear  relation  can  be  obtained  between  the  surface  and  >1  krn  by  asm: 
Eq.  (21)  and  holding  the  Junge  Slope  constant  at  -2  for  heights  a bov- 
4  km.  The  same  can  be  applied  to  the  summer  data  so  that  the  summer 
equation  becomes 


v  =  4  -  /  /  2  (2d) 

where  again  the  Junge  Slope  =  -  (  v  +  1  )  as  given  in  the  power  -  lav 

Eq.  (2),  and  Z  is  the  height  above  ground  in  km  and  can  not  exceeu  a 
value  of  4km.  At  heights  above  4km  in  summer  v,  equals  a  constant  value 
of  2  so  the  Junge  Slope  remains  constant  at  -3. 

We  now  want  to  obtain  an  equation  to  describe  the  decrease  in  ’he 
total  aerosol  concentration  with  height.  Fig.  9  shows  two  models  and  two 
data  sets  and,  although  they  vary  widely  m  the  value  for  particle  concen¬ 
tration  at  a  specified  level,  they  all  show  an  exponential  decrease  with 
height  from  the  surface  to  about  5km.  Above  that  height,  *otal  particle 
concentration  remains  essentially  constant  up  to  the  tropopause.  In  1954 
Penndorf  derived 

N2  =  N0  exp  (  -  Z  /  H  )  (23) 

to  give  the  total  number  concentration  of  aerosols  at  a  given  height,  Nz, 

as  a  function  of  the  total  concentration  at  ground  level,  N0,  the  height  in 
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the  atmosphere  in  km,  Z,  and  a  term  called  the  scale  height  in  km,  H. 
Penndorf's  data  suggested  variations  in  H  from  }  2  to  1.4  km.  More  re¬ 
cently  in  1970,  Elterman  made  simultaneous  observations  of  meteorological 
range  (  R^  ),  optical  attenuation  coefficients,  and  aerosol  scale  height.  Ir 
Fig.  13,  we  have  plotted  hts  data  points  and  derived 

H  =  0.8  +  Rm  /  30  (24) 

where  H  and  Rm  are  in  km.  This  Eq.  (24)  applies  for  visibilities  ranging 
from  0  to  18  km.  The  scale  height  remains  constant  at  1.4  km  for  visi¬ 
bilities  greater  than  18  km.  By  using  Equations  (23)  and  (24)  it  can  be 
seen  that  as  surface  visibility  increases  the  aerosol  concentration  with 
height  decreases. 

Thus  it  is  possible  to  use  the  relations  developed  in  this  section 
to  specify  number  density,  concentration,  and  ve  al  distribution  of 
aerosols  as  a  function  of  routinely  available  weather  data. 


Rm  =  METEOROLOGICAL  RANGE  1km) 


Pig.  13.  Aerosol  Scale  Height  Versus  Meteorological  Range 
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2. 1.1.3  Dynamic  Procedure 

Now  we  will  present  our  sequence  of  steps  taken  to  use  standard 
meteorological  observations  to  infer  continental  haze  particle  size  distribu¬ 
tions  having  radii  equal  to  or  greate.  than  0.1  micrometer.  These  steps 
are: 

1.  Determine  surface  air  particle  concentration  factor  C 
from  the  observed  meteorological  range  (Rm)  in  kilo¬ 
meters  where  C  -  30/  Rm. 

2.  Derive  the  surface  air  particle  size  distribution  using 

n  (Z  ,r)  =  =  0.434Cr"(v  +  *  )  where  n(Z ,r)  is  in  units 

of  cm'^  */m-1  and  v  is  an  exponential  factor  related  to 
the  type  of  aerosol  and  the  altitude  above  the  earth's 
surface.  We  found  the  haze  exponential  factor  as  a 
function  of  altitude  Z  in  kilometers  to  be  v  -  3  -  Z/2 
for  times  other  than  summer  months,  which  is  valid  from 
the  surface  to  4  km  where  v  remains  constant  a.  one  for 
higher  altitudes.  During  summer  months  v  =  4  -  Z/2 

1 

from  the  surface  to  4  km  above  which  v  remains  constant 
at  2 . 

3.  Determine  the  average  particle  concentration  N(o,r)  at 
the  earth's  surface  by  using  the  equations  for  C  and 
n(Z,r)  above  and  integrating  the  latter  for  spring  and  fall 
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conditions  from  r  to  infinity.  Thus 
N  (o,r)  =  0.  434  C  (l/3)r“3 
For  summer  conditions 

N  (o,r)  =  0.434  C  (l/4)r'4 

4.  Compute  the  decrease  in  average  particle  concentration 
N(Z,r)  as  altitude  (in  kilometers)  increases  above  th. 
ground  using  our  empirical  equation 

N(Z,r)  =  N(o,r)  exp  (  -  7/H) 

where  H  is  the  atmospheric  scale  height  in  km  that 

varies  with  meteorological  range  as  follows 

H  =  0.8  +  R  /3 0 
m 

where  H  reaches  a  constant  value  of  1.4  for  very  nigh 
visibilities  (i.e.  Rm  r  18  km). 

5.  After  deriving  the  average  particle  concentration  a:  a 
particular  height  N(Z,r),  the  size  distribution  value  n(7,r)  is 


n(z,r) 

=  v  r"  1 

N  (Z,r) 

or 

Z/2  )  r-2 

n(Z,r) 

=  (  3  - 

N  (Z.r) 

for  spring  and  fall 

and 

Z/2  )  r~ 2 

n(Z,r) 

=  (  4  - 

N  (Z.r) 

for  summer 

where  Z  in  the  right  hand  side  of  the  equations  can  never 
exceed  4  km. 

Following  the  above  steps,  we  made  computations  for  aerosol  size 
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distributions  at  1.8  and  6.1  km  to  correspond  with  Cress's  aircraft  ob¬ 
servations  in  western  Europe,  Fig.  14  and  15.  Cress  particle  size  mea¬ 
suring  equipment  was  restricted  to  detections  over  the  radius  range  from  0.2 
to  6  Jim.  Considerable  scatter  exists  in  his  observations  which  were 
taken  over  a  wide  variety  of  weather  conditions  ranging  from  clear  to 
overcast  skies  and  from  excellent  visibility  down  to  a  minimum  of  5  km. 

Now  assuming  the  surface  visibility  is  20  km  (  the  same  as  used  ir. 
the  AFGL  Rural  Model  ),  the  surface  air  particle  concemva  tion  fac:or  C 
equals  1.5  and  the  scale  height  equals  1.4.  So  the  average  particle 
concentration  N(0,r)  at  the  earth's  surface  is  for  spung  and  fall 
N(0,r)  =  0.434  (1.5)  ( ^-)  r~3 
N(0,0.1)  =  217  cm  3  for  particles  '  0.1  j/m 

N(0,0.2)  =  27.12  cm  3  for  20  km  visibility  and  f>artic!es  :  0.2  fjr:, 

N(0,0.4)  =  3.4  cm  3  for  particles  ?  0.4  j/m 
and  _3 

N(0,1.0)  =  0.22  cm  for  20  km  visibility  and  particles  ?  1.0  a n 
At  an  altitude  of  1.8  km  then 

N(1 . 8  ,r)  =  N(0,r)  e  ~  (tfr) 

N(l. 8,0.1)  60.0  cm-3  for  particles  z  0.1  j/n 

N(1.8(0.2)  7.50  cm-3  for  particles  -  0.2  j/m 

N(1.8,0.4)  0.94  cm-3  for  particles  ?  0.4  j jm 

N(l. 8,1.0)  0.06  cm  3  for  particles  >  1.0  fjm 
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then 

n(l .  8  ,r)  =  (  3  -  ~)  r"1  N(1.8,r) 

n(1.8,0.2)  =  78.8  cm-3  //m  *  for  r  =  0.2  /in 

n(l  .8, 1 .0)  =  0.13  cm-3  /im-*  for  r  =  1.0  fx rn 

At  an  altitude  of  6.1  km,  where  Cress  aircraft  observations  wore  taken 

over  western  Europe,  and  assuming  a  20  km  surface  visibility  then 

N(6.1,0.2)  =  0.35  cm~3  for  particles  t  0.2  /im 
N(6.  1  ,  1,0)  =  0.003  cm-3  for  particles  ;  1.0  /i m 

and 

n(6.  1,0.2)  =  1.75  cm-3  //m-*  for  r  =  0.2 
n(6. 1  ,  1.0)  =  0.003  cm-3  //m-1  for  r  -  1.0  /in 

Since  from  Eqs.  (1)  and  (2)  dN/d  log  r  =  (r/0.434)  n(Z,r)  then  at  1.-'  -.n 

dN/d  log  r  =  36.3  cm~3  //m“*  for  r  =  0.2  /im 

dN/d  log  r  =  0.3  cm-3  /im-1  for  r  -  1.0  /im 
then  at  6.1  km 

dN/d  log  r  =  0.81  cm~3  /im-*  for  r  =  0.2  //m 
dN/d  log  r  =  0.007  cm-3  /im-*  for  r  —  1.0  /im 
In  order  to  provide  a  means  for  comparison  with  our  GAC  Model,  plots 
are  shown  in  Figs.  14  and  15  for  Cress's  aircraft  data,  Blifford's  impactor 
data,  and  two  models  of  the  Air  Force  Geophysics  Laboratory  (APGL)  i.e., 
AFGL  tropospheric  Model  and  the  AFGL  Rural  Model  for  a  20  km  visibility. 

It  can  be  seen  that  the  GAC  Model  provides  the  best  fit. 


48 


tIN  d  log  r  (cm" 


1000 


10  00  h 


1  1  000  h 


0  100 


0  001 


|_  ALTITUOE  1  8  km 
GAG  Model 

20  km  visibility  +  \  + 

AFGL  RURAL  MODEL  ' 
VISIBILITY  '-20  km 

0010  —  -  AFGL  TROPOSPHERIC 

MODEL 

—  BLIFFORD  IMPACTOR  DATA  \ 

♦  AFGL  C  130  ROYCO  DATA  ^ 

J _ I _ 


RADIUS  (micrometers) 


Comparison  of  Aerosol  Models  at  6.1  km 


2.2  HYDROMETEORS 


2.2.1  Fog 

2. 2. 1.1  Ground  and  Upper  Level  Variations 

This  section  on  hydrometeors  will  cover  fog,  cloud,  and  rain  drop 
size  distributions.  An  overview  of  these  features  can  be  obtained  from 
Fig.  16  where  it  shows  about  a  12  order  of  magnitude  change  in  number 
density  and  about  half  that  change  for  representing  the  possible  range  m 
droplet  radii.  Often  in  the  real  world  there  is  a  number  of  possible  com¬ 
binations  of  hydro  -  and  litho-  meteors  occuring  simultaneously.  High 
pressure  regions  with  clear  skies  and  stable  boundary  layers  are  condusive 
to  good  radiational  cooling  and  fog  formation.  Also  it  is  not  unreasonable 
to  expect  trapping  of  haze  particles  in  the  stable  surface  air  layer  and 
have  a  mixture  of  both  haze  and  fog  particles.  Other  event  combinations 
occur,  such  as  drizzle  and  fog,  yet  we  tend  to  focus  our  modelling  efforts 
on  a  simple  event  representation.  That  is  also  done  here,  however,  pro¬ 
visions  are  made  and  results  shown  to  combine  outputs  from  multi  models 
to  better  simulate  all  reported  weather. 

Many  types  of  fog  exist  to  produce  large  spatial  and  temporal 
variations  in  "seeability.  "  Recent  improvements  and  technology  advances 
in  both  particle  sampling  and  sizing  instruments  and  observational  platforms 
needs  to  be  exploited  in  carefully  conceived  and  implemented  field  ex- 
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Thick  Fog 


Fig .  16.  Volume  Concentration  of  Water  Droplets  fay  Size 
(Counted  in  1  -ftm  Intervals  of  Droplet  Radius ) 


penmental  and  data  analysis  programs.  Such  data  are  basic  to  defining 
microphysical  structure  as  well  as  to  generating  synoptic  and  mesoscale 
models  as  done  here. 

Fog  conditions  often  go  through  three  distinct  stages,  i.e.  initial 
oscillatory,  steady  dense,  and  break  up  stage.  Particle  sizes  during  the 
initial  oscillatory  stage  usually  have  a  narrow  range  of  values,  the  most 
frequently  occuring  radius  is  small,  the  liquid  water  content  is  less,  and 
haze  and  fog  particles  combine  to  weight  the  size  distribution  to  smaller 
sizes.  Dense  persistent  fogs  are  older  in  character,  having  higher  liquid 
water  content,  broader  range  of  particle  sizes,  larger  mode  radius,  and 
greater  EO  attenuation.  The  break  up  stage  depends  greatly  upon  the  fog 
type,  formation,  density,  thickness,  and  physical  processes  producing  the 
break  up.  The  particle  size  distribution  of  breaking-up  fog  reverts  back 
more  to  a  haze  like  character  with  an  ever  decreasing  residual  peak  in 
the  concentration  of  the  larger  mode  radius.  Typical  fog  parameters  are 
shown  in  the  following  Table  2. 


Parameter 


Table  2.  Typical  Fog  Parameters 

Fog  type 

Radiation  (newl  Radiation  (old) 


mode  radius  (*fm)  4 

breath  of  size  distribution  narrowest 
vertical  depth  (m)  variable 


10 

broad 

200 


Valiev  Advection 
8  10 

broad  broadest 

150  400 
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The  time  period  that  usually  separates  new  from  old  fogs  is  on  the  order  of 
2  hours  after  formation.  Several  other  names  have  been  used  to  describe 
these  fog  types.  For  example  newly  formed  fogs  are  sometimes  called  "se¬ 
lective  fog."  Older  or  dissipating  fogs  may  be  called  "stable  fog"  or 
"evolving  fog."  Likewise  advection  fog  near  the  sea  shore  is  called 
"coastal  fog."  Often  "low  mountain  fog"  Is  called  valley  fog  and  stratus 
clouds  intersecting  a  mountain  are  called  "high  mountain  fog."  We  know 
the  preferred  mode  radius  of  fog  has  seasonal,  geographical,  and  meteo¬ 
rological  variations  that  need  better  definition.  Until  these  refined  da'a 
become  available,  for  our  purpose,  the  above  listed  mode  radius  will  be 
used. 

Liquid  water  content  of  a  fog  must  be  specified  in  order  to  define  the 
fog  density  or  concentration  of  droplets.  Sensors  do  exist  for  directly 
measuring  liquid  water  content  and  these  are  the  most  accurate  and  reliable 
but  are  not  too  feasible  within  a  battlefield  environment.  Horizontal  visi¬ 
bility  observations  are  available  from  routine  and  nonstandard  observational 
sources  and  may  be  used  to  infer  liquid  water  content  of  fog.  A  number  of 
experimental  studies  have  been  made  of  liquid  water  content  and  atmospheri 
visibility.  We  have  averaged  some  of  these  results  and  plotted  three  curve 
in  Fig.  17  to  depict  liquid  water  content  versus  visibility  for  coastal  fog 
and  two  inland  fog  types  (new  and  old).  The  equations  used  in  this  study 
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VISIBILITY  (STATUTE  MILES) 

Fig.  17.  Relationship  between  Visibility  and  Liquid  Water  Content  m  Fog 


to  derive  fog  liquid  water  content,  in  g  m.  ^  as  a  function  of  meteoro¬ 
logical  visual  range,  in  miles,  are 

Coastal  Fog  WL  -  (27.15  R^)'1,43  (25) 

Inland  Fog 

Old  WL  =  (41.7  Rm)"K5  (26) 

New  WL  =  (58.8  R^)  1‘“4  (37) 

By  comparing  information  from  Table  2  and  Fig.  17  it  can  be  seen 
that  the  narrowness  or  breath  of  drop  size  distributions  is  dependent  on  fog 
type  which  in  turn  is  related  to  liquid  water  content.  That  is  to  say  'ha- 
for  a  fixed  visibility,  liquid  water  content  is  higher  and  the  drop  size 
distribution  broader  for  coastal  than  for  either  inland  fog  types.  In  order  for 
modelling  techniques  to  be  representative  the  exponential  function  shape  fac'or 
a  must  vary  as  fog  density  varies  to  produce  a  size  distribution  which  is 
very  narrow  (  where  a  is  large  )  when  the  liquid  water  content  is  low  and 
very  broad  (  where  a  is  small  )  when  liquid  water  content  is  high.  We 
analyzed  a  wide  variety  of  droplet  size  distributions  to  obtain  the  following 
shape  factor  relationship 

a  =  1-1.4  In  WL  (28) 

where  is  the  fog  liquid  water  content  in  g  m-^.  For  computational 
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simplicity  and  so  that  our  equations  in  section  1  of  this  report  would  not 
have  to  increase  in  complexity,  we  placed  two  restrictions  on  the  value 
of  a,  first  it  can  never  be  less  than  one, and  second  it  must  be  rounded 
off  if  necessary  to  be  a  whole  integer. 

The  composition  of  inland  fog  normally  varies  with  height  above  ground 
level.  The  number  of  fog  droplets  per  unit  volume  usually  increases  from 
the  ground  to  the  top  of  the  fog.  Larger  droplets  and  higher  liquid  water 
content  are  usually  found  at  the  base  of  inland  fogs  whereas  coastal  fogs 
are  more  homogeneous  in  the  vertical.  For  older  and  more  stable  inland 
fogs  the  drop  size  distributions  become  more  narrow  and  unimodal  and  the 
mean  radius  decreases  with  increasing  altitude  above  the  surface.  For  a 
composite  of  old  inland  radiational  and  valley  fogs  we  found  averace  verti¬ 
cal  variations  for  mode  radius  rc  and  liquid  water  content  W'l  to  be 


d  rc  /  d  Z  =  -  1  jim  /  100  feet 

(29) 

and 

d  WL  /  d  Z  =  -  0.4  g  m-3  /  100  feet 

(30) 

with  a  minimum  value  of  4  jim  applicable  to  rc.  Observations  of  coastal, 
advection,  or  marine  fogs  often  show  nearly  constant  or  increasing  liquid 
water  content  with  increasing  height  in  the  fog.  For  this  paper  we  will 
assume  constant  conditions  prevail  from  the  surface  to  the  top  of  coastal 
fogs . 
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2.2. 1.2 


Dynamic  Procedure 

1.  Use  the  output  of  analyses  of  meteorological  surface  and  upper 
air  observations,  such  as  our  CFAS  (Cloud  Fog  Analysis  System) 
or  our  CIVAS  (Cloud/Icing/Vlsibility  Analysis  System),  to  Identi¬ 
fy  the  type,  age,  and  thickness  of  fog,  other  restrictions  to 
visibility  such  as  haze,  mist,  etc.  ,  and  the  visibility. 

2.  Select  the  applicable  mode  radii  rc  from  Table  2. 

3.  Derive  the  ground  level  liquid  water  content  from  either  Eqs. 

(25).  (26),  or  (27). 

4.  Compute  the  shape  factor  from  Eq.  (28). 

5.  Obtain  the  coefficient  A  from  Eq.  (19). 

6.  Derive  the  B  coefficient  from  Eq.  (10). 

7.  Determine  the  particle  size  distribution  n(r)  from  Eq.  (9). 

8.  Derive  the  total  nimber  concentration  of  droplets  for  all  radii 
using  Eq.  (l 1) . 

9.  Plot  the  output  of  steps  7  and  8  to  depict  surface  level  fog  drop¬ 
let  conditions. 

10.  Determine  if  other  constituents  are  also  restricting  the  visibility 
and  solve  for  and  incorporate  their  contributions  to  the  overall 
particle  size  distribution. 

11.  Derive  the  vertical  variation  of  rc  and  from  Eqs.  (29)  and  ( 3 J) 
and  incorporate  and  repeat  step  4  through  10  to  obtain  particle 
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size  distribution  at  any  desired  level  within  the  fog  layer. 

Using  the  above  steps  we  computed  particle  size  distribution  for 
observations  (Jiusto  1979)  of  haze  and  radiation  fog.  Meyer's  (1980)  da'.a 
showed  that  a  good  linear  relation  exists  between  visual  ranges  equal  to 
or  greater  than  5  km  and  the  cumulative  number  concentration  of  haze  parti 
cles.  The  number  concentration  almost  remains  constant  for  haze/fog  con¬ 
ditions  with  visual  ranges  from  1  to  5  km.  For  visual  ranges  less  than 
about  1  km  a  linear  relation  exists  between  the  number  concentration  and 
visual  range  in  dense  fog.  Therefore,  for  fog  and  haze  conditions,  the 
haze  component  of  particle  size  distribution  was  computed  for  a  visual 
range  of  5  km  and  was  held  constant  and  added  to  the  fog  contribution 
’•  ich  varied  with  visual  range.  Our  results  are  shown  in  Fig.  18  for  'wo 
reported  visibilities,  2.1  and  0,39  km.  For  comparison,  observed  size 
distributions  are  shown  for  three  observation  times  0725,  0730  and  0800 
and  visibilities  2.1,  1.4,  and  0.39  km,  respectively.  It  is  striking, 
however,  how  good  the  major  features  and  time  trends  are  represented  by 


our  haze/fog  model. 


0.05  0  5  5.0  50 


RADIUS  (micrometers) 


Fig.  18.  GAG  Haze-Foq  Model  Compared  With  Observed  Data 


Wi 


2.2.2  Cloud 


2.2.2. 1  Cloud  Base  and  In  Cloud  Variations 

We  have  had  a  great  deal  of  difficulty  looking  for  common  donomina'or; 
that  are  applicable  to  the  many  cloud  physics  studies.  This  is  due  in  par' 
to  the  fact  that  most  investigators  direct  their  attention  to  one  particular 
cloud  type  and  then  do  not  measure  all  variables  important  to  a  more  general¬ 
ized  study.  For  example.  Table  1  shows  that  for  a  particular  cloud  type  tr. are 
is  a  most  frequent  particle  radius  that  occurs  on  the  average  for  all  season-., 
and  types  of  conditions.  Yet  it  is  not  unreasonable  to  expect  that  cloud 
particles  and  their  distribution  are  dependent  upon  the  environment  in  which 
they  are  formed.  We  have  tried  to  take  these  features  into  accoun'  from  a 
meso-  and  synoptic-scale  point  of  view,  knowing  full  well  that  ccr'air.  micr- 
scale  features  must  either  be  neglected,  averaged,  or  inferred  from  currer.* 
observations . 

Our  desire  is  to  be  able  to  derive  cloud  characteristics  at  the  cloud 
base  and  heights  within  the  cloud,  given  only  standard  meteorological  ob¬ 
servations  . 

We  have  found  that  clouds  formed  primarily  by  convection  (cumulus 
and  cumulonimbus),  turbulence  (Stratus,  stratocumulus ,  and  altocumulus)  an.i 
horizontal  convergence  (altostratus ,  as  well  as  altocumulus)  exhibit  distinctive 
features  but,  most  often  when  looked  at  in  detail,  their  microphysics  is  dic- 
rated  by  the  temperature,  pressure,  dew  point,  and  vertical  motion  tha'  exists 
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at  the  cloud  base  and  levels  within  the  cloud.  That  is  to  say,  originally 
it  was  thought  that  each  cloud  type  would  have  to  be  trea'cd  separately. 

Now,  however,  we  found  that  cloud  base  characteristics  and  their  vertical 
variation  within  a  cloud  could  be  formulated  using  the  most  recent  surface 
and  upper  air  observation. 

Lewis  (1951)  compiled  tables  displaying  cloud  droplet  and  liquid  water 
content  for  a  large  number  of  aircraft  observations  segmented  m*o  *hree  hr  or. 
cloud  type  categories,  i.e.  stratus  and  stratocumulus ,  alto-stratus  and  alto¬ 
cumulus,  and  cumulus  or  cumulonimbus.  Separating  these  categories  geo¬ 
graphically  between  the  Pacific  Coast  and  other  regions  of  the  United  Stater 
showed  that  in  general  cloud  drops  for  corresponding  west  coast  cloud  types 
are  about  2  micrometers  larger  in  radius.  Although  cloud  base  temperature 
and  to  a  lesser  extent  the  cloud  base  pressure  height  contribute  to  chances 
in  the  mean  particle  size  we  were  unsure  whether  sufficient  differences  ex¬ 
isted  between  clouds  on  the  Pacific  Coast  versus  those  in  Eastern  U.5.A. 

We  used  data  prepared  by  the  Naval  Weather  Service  in  1976  on  "Climatic 
Study  of  the  Near  Coastal  Zone"  from  two  publications,  'East  Coast  of  the 
United  States"  and  "West  Coast  of  the  United  States."  On  the  average  the 
surface  air  temperature  is  10°F  warmer  in  the  winter  a'  San  Francisco  than  at 
Philidelphia  and  the  reverse  in  summer,  so  the  yearly  surface  air  temperatures 
are  about  the  same  value.  More  important,  however,  is  that  the  west  coast 
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clouds  occur  significantly  more  frequently  at  lower  heights  with  warmer  cloud 
base  temperatures.  This  is  especially  noticeable  during  the  winter  when  rr.os* 
of  the  clouds  exist  and  when  more  than  50  per  cent  of  the  west  coast  clouds 
have  a  ceiling  height  of  less  than  1000  feet  whereas  only  15  per  cent  are 
found  at  the  same  low  levels  on  the  east  coast.  We  also  analyzed  Selby's 
drop  size  measurements  in  low  level  stratus  in  another  country,  England, 
(Blifford  1970).  This  was  interesting  since  most  of  the  individual  cases  gave 
a  mode  radius  of  2  micrometers  for  cloud  droplets  within  the  lower  30  to  60 
meters  of  their  stratus  clouds.  Such  small  droplets  near  the  cloud  base  are- 
reported  more  frequently  in  the  later  literature  as  instrument  neasunng  tech¬ 
nology  improved.  The  main  point  is  that  major  features  of  these  liquid  clout. s 
from  wide  geographical  locations  were  found  to  be  represented  by  our  follow¬ 
ing  procedures.  Spatial  distribution  of  the  microstructure  of  cloud  liquid 
water  content  varies  considerably  in  both  horizontal  and  vertical  directions. 
This  is  especially  true  for  convective  type  storms  where  vertical  cloud  de¬ 
velopment  is  more  pronounced,  vertical  motions  are  higher,  the  total  cloud 
cellular  structure  is  often  composed  of  a  combination  of  individual  sub  cells, 
and  the  turbulent  motions  produce  more  entrainment  of  drier  ambient  air  which 
reduces  the  available  liquid  water  and  causes  gradients  in  the  actual  liquid 
water  content.  A  possibility  does  exist  to  use  remote  probing  techniques  to 
better  define  the  actual  cloud  microstructure.  However,  since  this  study  is 
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restricted  to  using  standard  surface  and  upper  au  observations,  we  r.ave  ■ 
develop  techniques  to  depict  major  features  associated  with  these  mc-s'.- 
scale  phenomena.  Using  information  contained  m  our  CIVAS  (cloud/' icing 
visibility  analysis  system)  we  can  specify  the  cloud  base  height  and  temper¬ 
ature,  cloud  type,  and  vertical  and  horizontal  extent.  From  this  we  car. 
compute  the  available  liquid  water  content  produced  by  cloud  air  rising 
along  its  moist  adiabatic  lapse  rate.  This  approach  produces  the  amount 
of  liquid  water  that  can  be  expected  at  each  level  in  the  cloud.  As  dis¬ 
cussed  above  some  of  this  water  has  to  go  to  injecting  moisture  into  the 
dry  entrained  air  to  bring  it  to  saturation  and  thus  the  actual  lie -:ti  water 
contained  at  any  cloud  level  is  less  than  that  expected  from,  purely  ac.a- 
batic  processes.  Furthermore  this  effect  and  moisture  reduction  ,ncrc 
with  altitude  in  the  cloud.  A  number  of  cloud  physics  studies  have 
made  showing  the  changes  of  the  ratio  of  actual  cloud  licuid  water  w  tha- 
expected  adiabatically  as  a  function  of  height  above  the  cloud  base, 
Pruppacher  (1980).  In  general  the  actual  liquid  water  is  a  high  petcentag-r 
of  that  available  near  the  cloud  base  and  decreases  to  nearly  a  quarter  of 
that  available  by  one  kilometer  above  the  cloud  base  and  remains  essentially 
constant  at  higher  cloud  levels.  We  have  fitted  these  data  with  a  linear 
curve  covering  the  first  kilometer  of  the  cloud  and  a  constant  value  of  C.2 
for  the  liquid  water  content  ratio  at  higher  levels  in  the  cloud.  The  equa¬ 
tions  are 
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WL  '  Wal  ~  1  ~  0,8  ZACB  for  ZACB  -  0  1  v-m  '31 

and 

WL  '  WAL  =  0-2  for  ?.ACB  >  ,  <rTi  (32 

where  Z  is  the  height  in  km  above  the  cloud  base  and  the  ra’m  is  *ne 
cloud  liquid  water  content  (Wr )  relative  to  the  cloud  liquid  water  con’ttv 
available  (WAjJ  through  moist  adiabatic  processes.  This  information  will 
be  combined  with  other  cloud  characteristics  to  obtain  a  measure  of  size- 
distribution  and  number  concentration  of  cloud  droplets  a-  different  heights 
in  the  cloud. 

Complicated  equations  are  necessary  to  precisely  derive  the  con:. licatf  d 
pseudoadiabatic  lapse  rate  as  a  function  of  temperature  and  pressure  at  tr .< 
cloud  condensation  level,  including  both  the  liquid  or  ice  stage  of  the  cloud. 
Another  set  of  equations  is  necessary  to  derive  the  saturation  mixing  ratio 
over  water  and  over  ice  as  a  function  of  air  temperature  and  pressure  with¬ 
in  the  clouds.  By  incrementally  solving  these  equations  it  is  possible  to 
derive  the  available  liquid  water  at  each  level  in  a  cloud.  This  precision 
is  unwarranted  at  this  time.  Considering  the  uncertainties  in  other  approx¬ 
imations  to  depicting  cloud  microphysical  features,  we  have  derived  the 
following  simp  ified  equations  to  derive  the  available  liquid  water  concentra¬ 
tion,  W^,  produced  by  moist  adiabatically  lifted  air 
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for  the  first  kilometer  in  the  cloud 


WAL  =  (  i  .42  *  0.05  ICB)  ZACB  ’33) 

and  for  heights  greater  than  one  kilometer  above  cloud  base 

WAL  =.  1.42  *  0.05  Tcb  *(0.84  -  0.035  ICE)  facr  1)  (34) 

a 

where  is  in  g  m  ,  ^C8  *s  t^'e  c^°uc^  base  temperature  in  °C  and  Z,^q  - 

is  the  height  above  cloud  base  in  km. 

Next  we  wanted  to  obtain  a  quantitative  method  to  derive  the  mod-'  n- 
dius  rc  since  this  is  a  necessary  variable  in  utilizing  the  previously  des¬ 
cribed  exponential  distribution.  The  peak  radius  of  a  size  distribution 
curve  was  found  to  be  directly  proportional  to  the  amount  of  liauid  water 
and  the  cloud  base  temperature.  We  empirically  derived 

r  =  1/5  (35) 

c  \340  -  8  tcb  / 

where  rc  is  the  mode  radius  in  micrometers,  is  the  cloud  liquid  water 
_  1 

content  in  g  m  ,  and  Tqb  is  the  cloud  base  temperature  in  °C .  This 
equation  was  derived  assuming  a  shape  factor  a  -  1  and  using  continental 
cumulus  cloud  data  to  derive  N  -  340  -  8  Teg  in  order  to  relate  the  aver¬ 
age  total  cloud  base  droplet  concentration  to  cloud  base  temperature.  Be¬ 
cause  of  the  limited  time  available  for  this  study,  this  equation  was  • non 
used  for  all  cloud  types.  In  order  to  derive  the  mode  radius  at  the  cloud 
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base,  we  assumed  no  entrainment  within  the  first  tenth  of  a  kilometer  of 
the  cloud  to  derive  W^,  which  is  then  dependent  only  upon  the  cloud  base 
temperature.  Therefore,  the  mode  radius  rc  at  any  cloud  base  is  given 
only  by  the  cloud  base  temperature,  producing  mode  radii  equal  to  2.8, 

3.5,  and  4.5  micrometers  for  cloud  base  temperatures  of  10,  20  and  30°C 
respectively. 

We  then  explored  two  methods  to  depict  vertical  changes  of  the  mode 
radius  within  a  cloud.  We  used  aircraft  observations,  Blifford  1970,  of 
summer  cumulus  and  obtained  the  following  best  fit  equation 

rc,Z  =  rc , Zcg  +  2‘3  ZACB  <36> 

where  rc<z  and  r c,Zq%  is  the  mode  radius  at  any  height  Z  within  the  cloud 
and  at  the  height  Z^g  of  the  cloud  base,  respectively,  and  Z^Qg  is  the 
height  in  kilometers  above  cloud  base.  The  other  method  was  to  use  the 
above  equations  to  derive  liquid  content  and  mode  radius  at  any  given 
height  within  a  cloud  of  known  base  temperature.  Also  used  was  Eq.  (23) 
to  derive  the  shape  factor  a,  Eqs.  (10)  and  (12)  to  derive  the  B  and  A 
coefficients,  and  Eq.  (11)  to  derive  the  total  droplet  concentration  N  per 
cubic  centimeter  volume.  Before  looking  at  detailed  comparisons  of  computed 
versus  observed  droplet  characteristics  as  a  function  of  height  within  a 
cloud,  we  will  look  at  comparisons  with  more  grossly  averaged  cloud  char¬ 
acteristics  . 
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Observations  from  5  different  investigators  (aufm  Kampe  i*  Weickmar.  HS7) 


were  combined  to  produce  a  frequency  distribution  of  mean  '.'.neai  i t " ■  r ! •" * 
radius  and  water  content  for  layer  (stratus)  clouds,  fair  wev.-.or  tumu'.us  , 
and  cumulus  congestus  or  cumulonimbus .  We  assumed  an  average  do-: 
base  temperature  of  15°C  and  used  the  observed  liquid  water  cornem  'o 


compute  the  expected 

mode  radius. 

The  res 

;ults  are: 

Observed 

Com  puted 

liquid 

mean 

mode 

mode 

difference 

Cloud  Tvoe 

water 

radius 

radius 

radius 

- 

stratus 

0.15 

5 .  S 

3.7 

2.8 

n  <  c 

fair  weather  cumulus 

0.6 

5.5 

3 .  7 

4.4 

0. 7 

cumulonimbus 

2.4 

12.5 

H.3 

6.9 

1.4 

The  internal  microphysical  structure  of  i 

layer  clouds  (aufm 

Ka  m  pe 

Weickman  1957)  was 

subjected  to  the  same 

analysis 

and  assumptions  as 

above  to  obtain  averaged  conditions 

for  the 

base,  middle,  and 

top  of  these 

clouds.  The  results 

are: 

Observed 

Computed 

Cloud  Type 

Location 

liquid 

mean 

mode 

mode  difference 

waiei 

radius 

radius  i 

radius 

stratus 

base 

0.03 

5.5 

3.7 

1 .  t>  2.1 

nimbostratus 

middle 

0.  15 

5.5 

3.7 

2.8  0.9 

top 

0.3 

8.4 

5.5 

3.5  2.0 

stratocumulus 

base 

0.15 

5.5 

3.7 

2.8  9.9 
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Observed  Computed 


Cloud  Type  (cont.) 

Location 

liquid 

water 

mean 

radius 

mode 

radius 

mode 

radius 

different.' 

altocumulus 

middle 

0.30 

5.5 

3.7 

3.5 

0.2 

(altostratus) 

top 

0.15 

5.5 

3.7 

CO 

M 

0.5 

It  is  interesting  that  the  differences  between  computed  and  observed  mode 
radius  are  about  one  micrometer  for  almost  all  cloud  types,  except  stratus 
where  computed  values  are  about  two  micrometers  too  small. 

We  looked  at  greater  vertical  detail  of  cumulus  cloud  observations 
(aufm  Kampe  &  Weickman  1957)  and  used  both  methods  described  above  to 
compute  expected  conditions.  In  this  case  the  average  cloud  base  temper¬ 
ature  was  25°C  and  we  obtained  the  following  comparisons  for  cumulus 
clouds: 


Height 

Within  Observed  Computed  "A"  Computed  "9" 


Cloud 

Lc 

Wt. 

N 

Ic 

Wt 

N 

r 

W^ 

N. 

base 

2.0 

0.  15 

330 

3.8 

0.25 

24b 

3.8 

0.25 

2-6 

1  km 

9.2 

1 . 0 

150 

4.8 

0.53 

153 

6.1 

0.53 

7  4 

3  km 

9.9 

2.8 

60 

6.4 

1.22 

50 

10.7 

1 .  22 

10 

5  km 

7.9 

2.0 

50 

7.4 

1.91 

52 

15.3 

1.31 

1  3 

where  computed  "B"  uses  empirical  Eq.  (36)  relating  mode  radius  as  a 
function  of  only  initial  conditions  at  cloud  base  and  height  above  it  where 
as  computed  "A"  uses  the  computed  vertical  distribution  of  liquid  water  and 
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the  corresponding  shape  factor  and  mode  radius  for  a  given  clouu  base  tem¬ 
perature.  In  all  cases,  total  particle  number  content  is  best  obtained  at 
all  levels  within  a  cumulus  cloud  by  using  method  "A",  whereas  in  nos- 
cases  method  "B“  provides  the  best  representation  of  the  mode  radius.  'he 
liquid  water  content  is  computed  the  same  for  both  methods  and  is  always 
somewhat  smaller  than  observed.  This  implies  our  entrainment  equation  is 
exerting  slightly  greater  influence  than  required  for  this  case.  Actually, 
with  a  higher  liquid  water  content  computed,  a  larqer  mode  radius  would  be 
computed  and  method  "A”  would  be  most  applicable  overall. 


2.2.2 .2  Dynamic  Procedure 


1.  Using  objective  surface  and  upper  air  analysis  techniques,  such 
as  CIVAS,  obtain  cloud  type,  cloud  base  temperature,  cloud 
height  above  ground  level,  and  horizontal  and  vertical  extent. 

2.  Compute  the  liquid  water  content  applicable  at  or  near  the  cloud 
base  using  Eq.  (33)  with  Zacb  =  0.1  km. 

3.  Combine  the  cloud  base  temperature  and  liquid  water  content  in 
Eq.  (35)  to  derive  expected  cloud  base  droplet  mode  radius. 

•4.  Derive  the  cloud  liquid  water  content  as  a  function  height  in  the 
cloud  by  applying  either  Eqs.  (33)  or  (34)  and  applying  -he  en¬ 
trainment  factor  Equations  (3l)  or  (32),  depending  upon  height 
above  cloud  base. 

5.  Determine  the  exponential  shape  factor  using  Eq.  (28)  at  each 
desired  height  in  the  cloud  where  liquid  water  content  was  com¬ 
puted. 

6.  Use  Eqs.  (9),  (10),  (11),  and  (12)  to  derive  the  details  of  the 
number  concentration  and  drop  size  distribution  as  a  function  of 
horizontal  and  vertical  distance  within  the  cloud. 
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2.2.3  Rain 


2.2.3. 1  Ground  and  Upper  Level  Variations 

Wartime  radar  research  provided  special  urgency  in  measuring  raindrop 

size  distribution.  In  1943  Laws  and  Parsons  began  a  new  approach  ‘o  ‘he 

problem  by  collecting  raindrop  size  data  and  relating  them  to  the  intensity 
of  precipitation.  They  found  as  rain  intensity 

increased  so  did  the  average  raindrop  size.  Also  they  could  use  climato¬ 
logical  rain  rate  data  to  infer  drop  sizes  and  effects  on  radar  for  differer.' 
geographical  locations  and  seasons  of  the  year.  Today  our  communications 
and  weather  reporting  network  is  such  that  daily  and  hourly  routine  meteoro¬ 
logical  observations  are  available  on  a  global  basis.  After  World  War  I! 
the  Stormy  Weather  Research  Group  at  McGill  University  studied  weather 
radar  responses  to  rainfall.  Marshall  and  Palmer  (1949)  then  found  that  the 
raindrop  size  distribution  relative  to  rain  rate  could  be  fitted  by 

N  (D)  =  Nq  exp  (  -  bD  )  (37) 

where  N  (D)  is  the  number  raindrop  per  unit  volume  (  m  )  and  per  unit  drop 
diameter  D  (mm),  NQ  is  the  limiting  value  of  N  (D)  for  D  =  0  and  is  often 
taken  as  NQ  =  8000  m-2  mm**  and  the  coefficient  b  in  units  mm-*  is  re¬ 
lated  to  the  rainfall  rate  (R)  in  mm  hr-*  by 

b  =  4.1  R-0-21  mm-1  (38) 
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Fig.  19  shows  the  raindrop  size  distribution  as  a  function  of  rainfall  rate, 


Laws  and  Parson  data  (broken  line),  and  observations  at  Ottawa  (dotted 
lines).  It  can  be  seen  that  good  correspondence  exists  for  all  raindror. 
sizes  and  rainfall  rates  except  for  the  smaller  drops  where  the  Marshall- 
Palmer  method  overpredicts  the  number  of  small  raindrops.  There  is  also 
an  upper  raindrop  size  limit  because  large  drops  (5  or  6  mm)  become  un¬ 
stable  and  break  up.  Cole  et  al  (1969)  suggests  limiting  the  Marshall- 
Palmer  method  to  raindrop  diameters  between  0.75  and  2.25  mm  for  rain 
rates  around  1  mm  hr"l,  between  1.25  and  3  mm  for  rain  rates  near  5  mm 
hr-*,  and  between  1.5  and  4.5  mm  for  rain  rates  greater  than  25  mm  hr-'. 
In  general,  however,  the  Marshall-Palmer  method  performs  reasonably  well 
to  provide  an  average  drop  size  spectrum  for  a  given  rainrate.  In  Switzer¬ 
land,  Joss  found  the  Marshall-Palmer  model  performed  well  for  the  same 
type  of  continuous  precipitation  in  which  it  was  developed  but  for  drizzle1 
and  thunderstorm  precipitation  the  coefficient  NQ  had  to  be  increased  anc: 
decreased  by  a  factor  of  4  and  l/8th,  respectively. 


Fig.  20  shows  the  large  variety  of  raindrop  size  distributions  that  exist 
for  different  goegraphicai  locations,  type  of  rain,  and  rainfall  intensity. 
Blanchard's  curves  1  -  3  are  for  Hawaiian  in-cloud  measurements  made  a'  or 
near  the  dissipating  edge  of  non-freezing  orographic  clouds,  while  carves 
4-7  represent  data  taken  at  the  cloud  base.  Curves  8-9  are  for  non- 
orographic  rain  distributions.  Curves  1-3  are  typical  of  what  is  expected 
in  a  combined  light  rain,  drizzle,  and  cloud  environment,  that  is  the  parti¬ 
cle  size  distribution  is  narrow  and  the  peak  frequency  occurs  at  very  small 
drop  sizes.  From  the  cloud  base  to  the  earth's  surface  it  can  be  seen  mat 
the  peak  frequency  or  mode  drop  diameter  increases  and  the  distribution 
broadens  as  rainfall  intensity  increases.  This  is  also  observed  in  our 
particle  size  and  number  concentration  distribution  for  mist,  drizzle,  light 
lain,  moderate  rain,  and  heavy  rain  in  Tables  i,  4,  5,  and  6,  respectively. 
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TABLE  3  Mist  and  Drizzle  Particles 


PARTICLE  SIZE 

WEATHER  TYPE 

MIST 

(0. 05  mm/hr) 


DRIZZLE 
(.25  mm/hr) 


AND  NUMBER  CONCENTRATION  DISTRIBUTION 

PARTICLE  NUMBER  PARTICLE  RADIUS 

_ DENSITY  (Km~3) _ (Km) 


6 

X 

109 

6.5 

X 

10"5 

13 

X 

109 

15 

X 

1  o"9 

51 

X 

109 

35 

X 

10"9 

69 

X 

109 

75 

X 

10"9 

51 

X 

109 

150 

X 

10"  0 

10 

X 

109 

250 

X 

1  0"  9 

2 

X 

109 

350 

X 

10"9 

.5 

X 

109 

450 

X 

ID"9 

.  1 

X 

109 

550 

X 

10"' 

6 

X 

109 

6.5 

X 

10"' 

13 

X 

109 

15 

X 

10"9 

51 

X 

109 

35 

X 

10"9 

80 

X 

109 

75 

X 

10"9 

85 

X 

109 

150 

X 

1  o~c_ 

26 

X 

109 

250 

X 

10  9 

9 

X 

109 

350 

X 

10 

3 

X 

109 

450 

X 

10"9 

1 

X 

109 

550 

X 

—  Q 

io  r 

.3 

X 

10 

650 

X 

10_9 

.  1 

X 

109 

750 

X 

10"9 

i 


,r 


CD  CD  O') 


TABLE  4  Light  Rain  Particles 


PARTICLE  SIZE  AND  NUMBER  CONCENTRATION  DISTRIBUTION 

PARTICLE  NUMBER  PARTICLE  RADIUS 

WEATHER  TYPE _ DENSITY  (Km~3) _ jKm) _ 


LIGHT  RAIN 


(1  mm/hr) 


6 

X 

I09 

6.5 

X 

10 

13 

X 

109 

15 

X 

10 

51 

X 

109 

35 

X 

10 

85 

X 

109 

75 

X 

10 

106 

X 

109 

150 

X 

10 

54 

X 

109 

250 

X 

10 

24 

X 

109 

350 

X 

10 

11 

X 

109 

450 

X 

10 

5 

X 

109 

550 

X 

10 

2 

X 

109 

650 

X 

10 

1 

X 

109 

750 

X 

10 

.4 

X 

109 

850 

X 

10 

.2 

X 

109 

950 

X 

10 

.3 

X 

109 

1250 

X 

10 
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TABLE  5  Moderate  Ram  Particles 


PARTICLE  SIZE  AND  NUMBER  CONCENTRATION  DISTRIBUTION 

PARTICLE  NUMBER  PARTICLE  RADIUS 

WEATHER  TYPE _ DENSITY  <Km~3) _ (Km} _ 


MODERATE  RAIN 


(4  mm/hr)  6  X  109 

13  X  109 
51  X  10® 
90  X  109 
146  X  109 
87  X  109 
49  X  109 
27  X  109 
15  X  109 
8  X  109 
4  X  109 
2  X  109 
1  X  109 
.2  X  109 
.1  X  109 


6.5 

X 

10 

15 

X 

10" 

35 

X 

10" 

75 

X 

10" 

150 

X 

10" 

250 

X 

10" 

350 

X 

1  0 

450 

X 

10" 

550 

X 

iC" 

650 

X 

10" 

750 

X 

10" 

850 

X 

10" 

9  50 

X 

10” 

1250 

X 

10 

17  50 

X 

10" 
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TABLE  6  Heavy  Rain  Particles 


PARTICLE  SIZE  AND  NUMBER  CONCENTRATION  DISTRIBUTION 

PARTICLE  NUMBER  PARTICLE  RADIUS 

WEATHER  TYPE _ DENSITY  (Km  3) _ (Kmj _ 


HEAVY  RAIN 


(16  mm/hr)  6  X 

13  X 
SI  X 
92  X 
160  X 
110  X 
75  X 
50  X 
32  X 
20  X 
13  X 
8  X 
6  X 
13  X 
2  X 
.1  X 


109 

6.5 

X 

10"9 

10® 

15 

X 

10~9 

109 

35 

X 

10'9 

109 

75 

X 

10"9 

109 

150 

X 

10"q 

109 

250 

X 

10'9 

109 

350 

X 

10_Q 

109 

450 

X 

io-9 

109 

550 

X 

10"9 

109 

650 

X 

10'9 

109 

750 

X 

IO-9 

109 

850 

X 

10"9 

109 

950 

X 

IO-9 

109 

1250 

X 

I  o'9 

109 

1750 

X 

IO-9 

109 

2250 

X 

10"9 
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We  have  developed  models  to  represent  the  simultaneous  vertical 


variations  in  rain  rate  and  liquid  cloud  content  from  the  earth's  surface 
to  cloud  top.  Shown  in  Table  7  and  8  is  our  summer  mid-latitude  mod¬ 
erate  (6  mm/ hr)  and  heavy  (15  mm/hr)  rain  model  which  provides  air  pres¬ 
sure  ,  temperature ,  relative  humidity,  cloud  content,  and  rain  rate  as  a 
function  of  height  in  the  atmosphere.  Notice  that  they  both  have  the  same 
cloud  base  and  top  but  differ  significantly  in  cloud  liquid  water  content  and 
rain  rate  outside  and  within  the  cloud.  The  rain  rate  decreases  rather 
slowly  (at  about  0.5  mm/hr/.  25  km)  from  the  surface  to  about  midway  through 
the  cloud  where  it  drops  to  near  zero  very  rapidly.  This  same  vertical 
variation  in  rain  rate  was  found  to  also  prevail  at  tropical  latitudes  for  the 
same  heavy  rain  situation,  Fig.  9,  but  where  the  cloud  top  extends  to  far 
greater  heights. 

The  general  form  of  the  exponential  function,  as  given  in  our  Eq.  (4), 
has  been  simplified  and  used  by  Bent,  Deirmendjian,  and  Khrgian  and 
Mazin  to  specify  static  models  of  rain.  If  we  take  our  Eq.  (17)  and  insert 
the  difference  in  velocity  between  the  updraft  (  Vu  )  and  terminal  fall  (  Vj  ) 
velocity,  (Vj  -  Vu),  in  the  integrand  it  transposes  the  left  hand  side  of 
the  equation  from  the  total  mass  concentration  of  liquid  water  into  the  rain¬ 
fall  rate.  When  the  updraft  is  small  relative  to  droplet  terminal  fall  velo¬ 
cities  then  the  equation  simplifies,  and,  as  shown  in  Fig.  21,  remains 


m — m - - ' 
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TABLE  7 


Moderate  Rain  (6  mm/hr)  For  Summer  Mid 


RAIN  MOOCl  ATNOS*K*t 


MOOCRAft  RAIN 


«  NN/M. 


rs  t  Moori  - 


MfICMT 

NtfUMI 

irM^fAATUHr 

PCLMIVt 

CLOUD  C04UNT 

IKNt 

(Mil 

toco  >) 

NUNIOITf 

(IN/CV  N» 

MM 

IRIS.RR 

IRI.» 

KIPP 

P.PPPO 

.IM 

Ral.RR 

2M.N 

l.PPP 

•  2PPP 

.Ml 

«A.RR 

TR5.IR 

I.PPP 

.2000 

•  tsp 

RIR.RR 

IRR.RR 

l.PPP 

•  32PP 

not 

HM.RR 

IR7.TR 

l.PPP 

.3300 

i.»o 

RTT.RR 

TRI.RR 

l.PPP 

.3400 

I.MI 

RAR.RR 

RAR.RR 

l.PPP 

•  35PP 

I.VSI 

RIR.RR 

ITR.SR 

l.PPP 

.3500 

MM 

m.ii 

ITT.1R 

I.POO 

•  35PP 

?»5M 

trt.rr 

ttr.rR 

l.PPP 

.3500 

MM 

TRI.RR 

tti.tr 

l.PPP 

.3500 

3*5P0 

RAR.RR 

7AR.5R 

l.PPP 

.3100 

OMI 

RIR.RR 

745 .RR 

l.PPP 

.2500 

4«MI 

STT.RR 

IRl.RR 

l.PPP 

.IPPO 

S«MI 

SrR.RR 

75R.RR 

l.PPP 

•  I3PP 

5.511 

SRS.RR 

7S4.5R 

I.PPO 

•  POPP 

6*000 

AW.RR 

751. SR 

I.PPO 

•  P5PP 

lilll 

RAR.RR 

IAA.7R 

l.PPP 

•  P3P0 

7. POO 

RIR.RR 

7RT.TR 

1.000 

.0100 

7.spp 

JRT.RR 

73R.RR 

1.000 

-*>.#000 

IrMI 

1AR.RR 

7TA.5R 

1  .000 

-P.PPPO 

o.stt 

11R.RR 

TTR.RR 

1.000 

-p.pppp 

I.MI 

IRT.RR 

TT5.RR 

.040 

-P.PPPO 

I.SM 

7<*5.RR 

TTR.RR 

•  050 

-p.pppp 

if. iff 

IRR.RR 

TTl.TR 

•  •50 

-•«•••• 

lUfff 

I’t.RR 

TIR.TR 

*4pp 

-p.pppp 

I2.PPP 

IRl.RR 

TIR.TR 

•  160 

-•••POP 

13*000 

1AS.RR 

TIR.TR 

•  POP 

-•••POP 

1 4*000 

I.I.RR 

TIR.TR 

.041 

-P.PPPO 

15. POP 

ITR.RR 

TIR.TR 

.031 

•p.pppp 

lf*Mf 

IRl.RR 

TIR.TR 

.•23 

-•.•••• 

17.111 

•  T.RR 

TIR.TR 

•  PIP 

-p.pppp 

11*111 

75. RR 

TIR.TR 

.016 

-p.pppp 

10.000 

RR.IR 

TIR.TR 

•  014 

-p.pppp 

20*000 

5R.TR 

TIR.TR 

.012 

-••PPPP 

12.111 

RR.R* 

TIR.TR 

.007 

-•*•••• 

MrMI 

ta. 

TIR.TR 

•  004 

-•••POP 

24.000 

71  .5R 

TTT.TR 

.002 

-p.pppp 

21 *11 P 

15.  *R 

TIR.TR 

.001 

-P.PPPO 

30*000 

II.  TR 

TIR.TR 

-0.000 

-p.pppp 

35*000 

5.RR 

717. RR 

-••POP 

-p.pppp 

4P.PPP 

T.RR 

7AI.RR 

-p.ppp 

-p.pppp 

45.PPP 

I.AR 

7R5.RR 

-p.ppo 

-p.pppp 

SP.PPP 

.RR 

TTR.TR 

•p.ppp 

-P.PPPO 

-Latitudes 


NlO-tATITUOf 

■AIM  RA»f 

(MM /MR  I 


R.RRR 

5.SRR 

A.AM 

).MI 

(.AM 
I.MI 
I.MI 
1  •  1  •• 

.At* 

-R.RRR 

-R.RRR 

-R.RRR 

•R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

•R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

•R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 

-R.RRR 
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TABLE  3  Heavy  Rain  (15  mm/hr)  For  Mid-Latitudes 


TABLE  9  Heavv  Rain  (15  mm/hr)  For  Troolcal  Latitudes 


UIOPICAL  -  15  M 

w 


.««« 


11. 


HC»»T  MIN  Of  19  HH/rtft 


mi w  m 


r.M 

i: ::: 

1.000 

l.ooi 


0#» 
i.ooo 

hoot 

i:s°oi 

i’.SoS 

i:IS5 
1.00  0 
1.00  0 
1.000 

i:coo 
i.cuo 
1.00  0 
I.o«e 

1.  Ov  0 

>00  0 


l.  10  0 
•  00" 


00  0 

it;! 

•  0<fi 

•  too 

:BI 

-0.080 


:!:!!! 


up  CONTCMT 

CU  hi 


•■•Mil 

ami 


[VM\ 

19  ' 


iS  :!l 


ViilU 
12.200 
11. too 
1 1.000 

»!  ' 


11.000 

10.900 


10.200 
9.90  0 
9.200 

VAU 

.200 


.  100 
-0.00* 


..  >00 

•  o.ooo 
-o. ooo 

-0.000 

-o.ooo 

-0. OkO 

:I:IU 

A:  m 
-8.000 
-0.  too 


:1V-'- 


.000 
-0.000 
-0.000 
-0.000 
-1.000 
-0.000 


-0.000 

-0.000 


-0.000 


-0.0000 


ifclH 


-0.000 


only  dependent  upon  raindrop  size.  This  rain  rate  equation  could  be  used 
as  a  feedback  mechanism  and  interate  between  it  and  the  assumed  form  of 
the  n(r)  equation.  That  is  to  say,  for  the  Marshall-Palmer  equation  one 
could  assume  an  initial  value  for  N  ,  solve  the  integral  to  obtain  a  first 
guess  value  for  R,  compare  the  first  R  value  with  that  observed,  select  a 
new  value  for  NQ  to  reduce  the  difference  between  observed  and  computed 
R,  and  iterate  to  obtain  the  most  applicable  N0  for  that  rain.  A  similar 
approach  can  be  used  with  the  exponential  or  modified  gamma  function  as  an 
aid  in  modifying  one  of  the  parameters  such  as  the  shape  factor  and/or  the 
coefficient  A  which  is  related  to  the  liquid  water  content.  Fig.  22  shows  th< 
log  log  relation  that  exists  between  rainfall  rate  R  (mm  hr-’)  and  liquid  w  - 
ter  content  {g  m-^)  of  raindrop.  Two  equations  express  this  analy'ica.ly 
as 

WL  =  0.072  R*88  (  ') 

for  Marshall-Palmer  rain  and 

WL  =  0.  1 1  R’88  (.  (■) 

for  joss  drizzle. 

In  fact,  for  many  obscuration  problems  it  is  possible  to  infer  sensor  atter.u- 

I 

ation  directly  by  meteorological  measurement  of  rainfall  rate  on  liquid  water 
content.  This  is  especially  noticeable  in  the  microwave  portion  of  the  spec¬ 
trum  as  shown  in  Fig.  23  as  well  as  in  the  visible  and  infrared  (Low  1979) 
for  other  weather  features  such  as  fog.  Here,  however,  we  are  interested 
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DROP  DIAMETER  (CM) 


Fig.  21  Terminal  Velocity  Vs.  Raindrop  Size 


LIQUID  WATER 


RAINFALL  RATE  (mmh'1) 


Fig .  2  2  Liquid  Water  Content  of  Precipitation  Versus  Rainfall  Rate 
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in  specification  of  raindrop  size  distributions  and  we  have  two  strong  can¬ 
didates,  the  exponential  approach  that  we  previously  used  for  fogs  and 
clouds,  and  the  Marshall-Palmer  method.  During  the  course  of  analyzing 
both  approaches,  we  analyzed  Marshall-Palmer  type  data  from  several  sour¬ 
ces,  geographical  locations,  and  rainfall  types  and  was  able  to  obtain  a 
best  fit  equation  that  solves  the  coefficient  N0  problem  previously  discussed. 
Figure  24  shows  our  plot  of  the  variation  of  NQ  coefficient  relative  to  rain¬ 
fall  rate  with  best  fit  equation  being 

Nq  =  12000  R"*72  (-11) 

This  equation  alters  the  raindrop  size  distribution  properly  by  increasing 
the  population  of  small  droplets  in  drizzle  while  decreasing  the  population 
of  raindrops  associated  with  showers  and  thunderstorms.  By  combining 
Equations  (37),  (38),  and  (41)  it  is  now  possible  to  have  a  single  analytical 
expression  for  deriving  raindrop  particle  size  distributions  for  all  rain/ 
drizzle/mist  type  environments  and  geographical  locations  using  only  routinely 
measured  rate  of  rainfall. 
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Fig.  24  Coefficient  N  Versus  Rainfall  Rate 
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Dynamic  Procedure 

1.  Using  objective  surface  and  upper  air  analysis  techniques,  suit 
as  CIVAS,  obtain  rainfall  rate,  height  of  cloud  base  and  top, 
and  horizontal  changes  in  precipitation  and  cloud  features. 

2.  Solve  Equations  (3  8)  and  (41)  to  obtain  applicable  values  of  : 
and  N0  given  the  observed  rainfall  ra*e  at  the  earth's  surface. 

3.  Solve  Eq.  (37)  to  obtain  raindrop  size  distributions  over  a  seiotr- 
ed  range  of  raindrop  diameters  from  0.7  5  to  2.2  5  mm  for  rainram:.- 
equal  to  or  less  than  1  mm/hr,  from  1.25  to  3  mm  for  rainrates 
greater  than  1  mm/hr  but  less  than  25  mm/hr,  and  from  1.5  to 
4.5  mm  for  rainrates  equal  to  or  greater  than  25  mm/hr. 

4.  Obtain  the  vertical  distribution  of  rainfall  rate  by  decreasing  the 
surface  value  by  0.5  mm/hr  per  0.25  km  height  in  the  a tmaspne:  <' 
up  to  the  mid-point  height  in  the  cloud  at  which  the  ramrato  .imps 
to  zero. 

5.  Use  the  above  derived  rainrate  distribution  in  the  horizontal  and 
vertical  to  derive  changes  in  raindrop  particle  size  distributions 
in  the  cloud-free  air  and  cloud  environment. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


1.  Dynamic  models  are  presented  for  specifying  par  le  size  dis¬ 
tributions  in  horizontal  and  vertical  directions  for  haze,  fog, 
cloud  and  rain  conditions  given  only  routinely  available  meteoro¬ 
logical  data.  Case  studues  showed  these  models  performed  well 
individually  and  when  used  together.  This  effort  should  be 
viewed  as  the  first  step  in  developing  dynamic  models  that  are 
responsive  to  observed  and  forecast  weather  changes.  These 
models  need  to  be  applied,  tested,  modified,  and  improved. 

2.  More  complete  microphysical  observations  are  needed  not  only 
to  better  understand  atmospheric  processes  but  also  to  provide 
better  inputs  to  such  meso-  and  synoptic  scale  models  as  de¬ 
veloped  here.  Efforts  should  be  made  to  incorporate  and  combine 
these  dynamic  models  with  such  Army  Cloud  Fog  Analysis  System 
(CFAS)  and  Cloud/Ice/Visibility  Analysis  System  (CIVAS)  to  de¬ 
pict  natural  obscurants  at  any  desired  time  or  location. 

3.  A  more  cooperative  working  environment  must  be  created  between 
the  micro-  and  macro-  atmospheric  physicist.  The  micro-physic  is 
feels  threatened  to  think  that  Army  users  could  be  satisfied  with 
only  macro-  scale  data.  In  turn,  the  macro-  physicis*  feels 
threatened  when  told  only  extensive  microphysical  observations 
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can  provide  necessary  details.  We  have  a  lot  to  accornphs: 
we  are  to  someday  be  in  a  position  to  help  our  Army  pr  jvid< 
day  to  day  and  hour  by  hour  assessments  of  defensive  and  of 
fens ive  weapon  effectiveness.  This  requires  all  talents  opera", 
collectively  with  realistic  guidelines  on  what  meteorological  d 
will  be  available  for  use  in  specifying  and  predicting  obscure" 
for  a  Field  Army. 


92 


|r  AD-A107  637 
UNCLASSIFIED 


GEO-ATMOSPHERICS  CORP  LINCOLN  Ma  p/G  20/6 

MODELS  OF  WEATHER  ENVIRONMENTS  ADVERSE  TO  ELECTRO-OPTICAL  SYSTE— ETCCU) 
SEP  81  W  D  MOUNT#  B  R  FOW  DAAD07-80-C-0047 

ERADCOM/ASL-CR-81-0047-1  NL 


4 


REFERENCES 


Berry,  E.W,,  and  R.L.  Reinhardt,  "Modeling  of  Condensation  and 
Collection  Within  Clouds,"  Phys ,  Sci.  ,  H>,  (1973). 

Blanchard,  D.C.,  and  Spencer,  A.T.,  "Experiments  on  the  Generation 
of  Raindrop-Size  Distributions  by  Drop  Breakup,"  J.  Atmos.  Sci.,  27 . 
(1970). 

Blifford,  I.H.,  Jr.,  et  al,  "Particulate  Models:  Their  Validity  and 
Application,"  NCAR  Technical  Notes,  Cat.  No.  NCAR-TN/PROC-68 ,  Nat. 
Center  for  Atmos.  Res.,  Boulder,  CO,  (1970). 

Blifford,  I.H.,  and  Ringer,  L.D.,  "The  Size  and  Number  Distribution  of 
Aerosols  in  the  Continental  Troposphere,"  J.  Atmos.  Sci.,  26.,  /\  (  j%9). 

Blifford,  I.H.,  Jr.,  "Tropospheric  Aerosols,"  J.  Geophys.  Research,  75., 

15,  (1970). 

Berry,  E.X.  ,  "Cloud  Droplet  Growth  by  Collection, "  J.  Atmos.  Sci.  24 . 

6,  (1967). 

Best,  A.C.  "The  Size  Distribution  of  Raindrops  ,"  Q.J.  R.  M.S.  ,  _76,  (1950). 

Bigg,  E.K.,  "Comparison  of  Aerosol  at  Four  Baseline  Atmospheric  Monitoring 
Stations,"  J.  Appl.  Meteorol.  ,  J_9,  (1980), 

Bigg,  E.K.,  "Some  Properties  of  the  Aerosol  at  Mauna  Loa  Observatory," 

J.  Appl.  Meteorol.,  Jji,  (1977). 

Brown,  R.  and  Roach  W.T.  "The  physics  of  Radiation  Fog:  II  -  A  Numerical 
Study,"  Q.J.R.M.S.,  102  ,  432,  p.  355,  (1976). 

Brown,  R.  ,  "A  Numerical  Study  of  Radiation  Fog  with  an  Explicit  Formulation 
of  the  Microphysics,"  Q.J. R. M.S. ,  1 06 .  450,  (1980). 

Carbone,  R.E.  and  Nelson,  L.D.,  "The  Evolution  of  Raindrop  Spectra  in 
Warm-Based  Convective  Storms  as  Observed  and  Numerically  Modeled ," 

15,  (1978  ) . 

Carrier,  L.W.  et  al.  "The  Backscattenng  and  extinction  of  Visible  and 
Infrared  Radiation  by  Selected  Major  Cloud  Models,"  Appl.  Opt.,  6., 

(1967). 


93 


Carswell,  A.I.,  "Laser  Measurements  in  Clouds,  "  IFAORS  Abstrac*  i>igcs- , 
Clouds:  Their  Formation,  Optical  Properties  and  Effects,  (1980). 

Chagnon,  C.W.,  and  Junge,  C.E.,  "The  Vertical  Distribution  of  Sub- N‘:  .;io 
Particles  in  the  Stratosphere."  Handbook  of  Geophysics  and  Space  Environ¬ 
ments,  J.  Meteorol.  JJ3,  p.  746,  (1961). 

Chappell,  C.F.  ,  and  Smith,  D.R.,  "Simulation  of  Cold  Cloud  Precipitation 
in  a  Three-dimentional  Mesoscale  Model,”  NOAA  Tech.  Rept.,  ERL-361- 
APCL-40 ,  26  pp.  ,  (1976). 

Chen,  C.C.,  "Attenuation  of  Electromagnetic  Radiation  by  Haze,  Fog, 
Clouds  and  Rain,"  USAF  Project  Rand,  R-1964-PR,  Santa  Monica,  CA.  , 
(1975). 

Christian,  H.,  Holmes,  C.R.,  Bullock,  J.W.,  Gaskell,  W.  ,  lllingsworti. , 
A. I.,  and  Latham,  J.  ,  "Airborne  and  Ground-Based  Studies  of  Thunderstorm, 
in  the  Vicinity  of  Langmuir  Laboratory ,  "  Q.J.R.M.S.,  1  06  .  (1980). 

Clark,  W.E.,  and  Whitby,  K.T.,  "Concentration  and  Size  Distribution 
Measurments  of  Atmospheric  Aerosols  and  a  Test  of  the  Theory  of  Self- 
Preserving  Size  Distributions,"  J.  Atmos.  Sci.,  2 4  ,  (1967). 

Cohen,  I.D.,  Capt,  USAF,  "Marine  Boundary  Layer  Sampling  :  Fliqh;  No. 

2,"  AF  Geoph.  Lab.,  AFGI.-TR-79-0242 ,  Res.  Papers,  678,  (1979). 

Cole,  A.E.,  Donaldson,  R.J.,  Dyer,  R.  ,  Kantor,  A.J.,  Skrivanek, 
"Precipitation  and  Clouds:  A  Revision  of  Chapter  5,  Handbook  of  Geoph ysic 
and  Space  Environments,"  USAF  Off.  of  Aerospace  Res.,  AFCRL-69-0497  , 
(1969). 

Cooper,  W.A.  ,  and  Saunders,  C.P.R.,  "Winter  Storms  over  the  San  Juar. 
Mountains,  Part  II:  Microphysical  Processes,"  J.  Appl.  Meteor.,  19 . 

(1980). 

Corradini,  C.,  and  l'onna ,  G.,  "The  Parameterization  of  the  Gravitational 
Water  Flux  in  Fog  Models,"  J.  Atmos.  Sci.,  .37  ,  (1980). 

Cotton,  W.R.,  “Cloud  Physics:  A  Review  for  1975-1978  IUGG  Quadrennial 
Report,"  Reviws  of  Geophysics  and  Space  Physics,  1 7  .  7  ,  (1979). 


94 


Cress,  T.S.,  "Altitudinal  Variations  of  Aerosol  Size  Distributions  Over 
Northern  Europe,"  AF  Geophys.  Lab.  Report,  Hanscom  AFB,  MA,  (1979). 

•.since  published  as  AFGL-TR-80-0178) 

Deirmendjian,  D.  ,  "Complete  Microwave  Scattering  and  Extinction  Properties 
of  Polydispersed  Cloud  and  Rain  Elements,"  USAF  Project  Rand  Report, 
R-422-PR,  (1963). 

Delany ,  A.C.,  Pollock,  W.H.,  and  Shedlovsky,  J.  P. ,  "Tropospheric  Aerosol: 
The  Relative  Contribution  of  Marine  and  Continental  Components,"  J.  Geophys. 
Res.,  23..  27,  p.  6249-6265,  (1973). 

Donnadieu,  G.,  "Comparison  of  Results  Obtained  with  the  VIDIhZ  Spectro- 
Pluviometer  and  the  Joss-Valdvogel  Rainfall  Disdrometer  in  a  "Rain  of  a 
Thundery  Type","  J.  Appl.  Meteor.,  J_9,  (1980). 

Downs,  A.R.,  "A  Review  of  Atmospheric  Transmission  Information  in  the 
Optical  and  Microwave  Spectral  Regions,"  USA  Ballistic  Res.  Labs, 

Aberdeen  Proving  Ground,  MD,  AD/A-035  059,  (1976). 

East,  T.W.R.,  "An  Inherent  Precipitation  Mechanism  in  Cumulus  Clouds," 
Q.J.R.M.S.,  83  ,  (1957). 

Eldridge,  R.G.,  "Haze  and  Fog  Aerosol  Distributions,"  J.  Atmos.  Sci., 

23  ,  (1966). 

Elterman,  L.  ,  "Parameters  for  Attenuation  in  the  Atmospheric  Windows  for 
Fifteen  Wavelengths,"  Appl.  Opt.,  _3«  745-749,  (1964). 

Elterman,  L.  ,  "UV,  Visible  and  IR  Attenuation  for  Altitudes  to  50  Km,  1968," 
AF  CRL-68-01 53 . 

Elterman,  L.  ,  "Vertical-Attenuation  Model  With  Eight  Surface  Meteorological 
Ranges  2  to  13  Kilometers,"  Hanscom  AFB  Report,  (1970). 

Falcone,  V.J.,Jr.,  Abreu,  L.W,,  and  Shettle,  E.P.,  "Atmospheric  Attenuation 
of  Millimeter  and  Submillimeter  Waves:  Models  and  Computer  Code,"  Hanscom 
AFB,  MA,  Report,  AFGi-TR-79-0253  ,  Optical  Phys .  Div.,  Proj .  7670,  (1979). 

Federer,  B.  ,  and  Waldvogel,  A.,  "Hail,  and  Raindrop  Size  Distributions 
from  a  Swiss  Multicell  Storm,"  J.  Appl.  Meteor.,  14,  p.  91-97,  (1975). 

Fenn ,  R.W.,  and  Shettle,  E.P.  ,  "Optical  Aerosol  Models  and  Light  Scattering 
Programs,"  AF  Geophys.  Lab.,  Hanscom  AFB,  MA . 


95 


Fenn,  R.W.,  Gerber,  W.  E. ,  and  Wasshausen,  E.(  "Measurements  of  the 
Sulfate  and  Ammonium  Components  of  the  Artie  Aerosol  of  the  Greenland  Ice 
Cap,'*  Handbook  of  Geophysics  and  Space  Environments,  (1963). 

Fitzgerald,  J.W.,  "Approximation  Formulas  for  The  Equilibrium  Size  of  an 
Aerosol  Particle  as  a  Function  of  Its  Dry  Size  and  Composition  and  the 
Ambient  Relative  Humidity,"  J.  Appl.  Meteor.,  .14.,  ( 1 975 ) . 

Fitzgerald,  J.W.,  "Angular  Truncation  Error  of  the  Integrating  Nephelometer 
in  the  Fog  Droplet  Size  Range,"  J.  Appl.  Meteor.,  1 6 ,  (l°77). 

Fletcher,  N.H.,  The  Physics  of  Ralnclouds,  Cambridge  University  Press, 
London,  pp.  112-121,  267-268,  (1966). 

Fujiwara,  M.  ,  "Raindrop-Size  Distribution  from  Individual  Storms,  "  J.  At:. 
Sci .  ,  22 ,  (1965). 

Garland,  I. A.,  "Some  Fog  Droplet  Size  Distributions  Obtained  by  an 
Impaction  Method,"  Q.I.R.M.S.,  92,  pp.  483-494  ,  (197  1). 

Gerber,  H.E.,  "A  Saturation  Hygrometer  for  the  Measurement  of  Relative 
Humidity  Between  95  and  105%,"  J.  Appl.  Meteor.,  19_,  (1980). 

Gillespie,  J.R.,  and  List,  R.  ,  "Evolution  of  Raindrop  Size  Distributions 
in  Steady  State  Rainshafts,"  Dept,  of  Phys.  ,  Univ.  of  Toronto,  (1976). 

Griggs,  M.  and  Marggraf,  W.A.,  "Measurment  of  Cloud  Reflectance 
Properties  and  the  Atmospheric  Attenuation  of  Solar  and  Infrared  Energy," 
Office  of  Aerospace  Res.,  Bedford,  MA,  Contract  No.  AF1 9  (628)- 5  5 1  7  , 
(1967). 

Hamilton,  P.M.,  "Vertical  Profiles  of  Total  Precipitations  in  Showtr 
Situations,"  Q.J.R.M.S.,  9 2_,  (1966). 

Hanel,  G.,  "The  Single-Scattering  Albedo  of  Atmospheric  Aerosol  Particles 
as  a  Function  of  Relative  Humidity,"  J.  Atmos.  Sci.,  3.3  ,  (1976). 

Herzegh,  P.H.,  and  Hobbs,  P.V.,  "The  Mesoscale  and  Microscale 
Structure  and  Organization  of  Clouds  and  Precipitation  in  Midlatitude 
Cyclones.  II:  Warm-Frontal  Clouds,"  J.  Atmos.  Sci.,  32,  3,  (1980). 

Hildebrand,  P.H.,  "Iterative  Correction  for  Attenuation  of  5  cm  Radar  in 
Rain,"  J.  Appl.  Meteor.,  \J_,  (1978). 


Hobbs,  P.V.  ,  Politovich,  M.K.,  and  Radke,  L.F.,  “The  Structures  of 
Summer  Convective  Clouds  in  Eastern  Montana,  I:  Natural  Clouds,"  J. 

Appl.  Meteor.,  19 ,  (1980). 

Hobbs,  P.V.,  Matejka,  T.J.,  Herzegh,  P.H.,  Locatelli,  J.D.,  and  Houze, 
R.A.  ,  ,rr.  ,  11  The  Mesoscale  and  Microscale  Structure  and  Organization  of 
Clouds  and  Precipitation  in  Midlatitude  Cyclones.  I:  A  Case  Study  of  a 
Cold  Front,"  J.  Atmos.  Sci.  ,  37.,  3,  (1980). 

Kobbs,  P.V.,  "Scales  Involved  in  the  Formation  and  Organization  of 
Clouds  and  Precipitation,"  IFAORS  Abstract  Digest,  Clouds:  Their  Forma¬ 
tion,  Optical  Properties  and  Effects,  (1980). 

Hofmann,  D.J.  et  al.,  "Stratospheric  Aerosol  Measurements  IV:  Global  Time 
Variations  of  the  Aerosol  Burden  and  Source  Considerations,"  J.  Atmos.  Sci., 
33,  9,  p.  1782-1788.  (1976). 

Hofmann,  D.J.  and  Rosen,  J.M.,  "Balloon  Observations  of  the  Time  Devel¬ 
opment  of  the  Stratospheric  Aerosol  Event  of  1974-1975,"  J.  Geophys.  Res., 
82.,  9,  p.  1435-1440,  (1977). 

Houze,  R.A.  ,  Hobbs,  P.V.,  Herzegh,  P.H.,  and  Parsons ,  D.B.,  “Size 
Distribution  of  Precipitation  Particles  in  Frontal  Clouds,"  J.  Atmos.  Sci., 

36  ,  (1979). 

Howell,  W.E.,  "The  Growth  of  Cloud  Drops  in  Uniformly  Cooled  Air," 

J.  Meteor.  ,  6.,  (1949). 

Hudson,  J.G.,  and  Squires,  P. ,  "Microphysics  of  Post  Ana  Fogs," 
International  Conference  on  Cloud  Physics,  Boulder  Colorado,  (1976). 

Hudson,  J.G.,  "Relationship  Between  Fog  Condensation  Nuclei  and  Fog 
Microstructure,"  J.  Atmcs.  Sci.,  2Z.«  (1980). 

Hutfaker,  R.M.,  et  al.,  "Feasibility  Study  of  Sattelite-Borne  Lidar  Global 
Wind  Monitoring  System,"  NOAA  Technical  Memorandum  ERL  WPL-37,  (1978]. 

Hunt,  G.,  "Diurnal  Properties  Of  Clouds  From  Geostationary  Satellite 
Observations,"  IFAORS  Abstract  Digest,  Clouds:  Their  Formation,  Optical 
Properties  and  Effects,  (1980). 


97 


[ 

I 

I 


Jones,  D.M.A.,  "Rainfall  Drop  Size  Distribution  and  Radar  Reflectivity,  " 
Meteor.  Lab.  Res.  Rep.  No.  6,  Illinois  State  Water  Survey  (NTIS  No. 

API  01  799.),  (1956). 

Johnson,  D.B.,  "The  Influence  of  Cloud-Base  Temperature  and  Pr«-ssur< 
on  Droplet  Concentration,"  J.  Atmos.  Sci.  3~ ,  (1980). 

Joss,  J.  ,  and  Waldvogel,  A.,  "Raindrop  Size  Dis'ribution  and  Samplin’: 

Size  Errors,"  J.  Atmos.  Sci.  ,  .26,  (l%9). 

Joss,  J.  ,  and  Gori,  E.G.,  "Shapes  of  Raindrop  Size  Distributions ,  " 

Appl .  Meteor.,  17  .  (1978). 

Jones,  M.A.  ,  and  Sims,  A.  L. ,  "Climatology  of  Instantaneous  Rainfall 
Rates,"  J.  Atmos.  Sci.,  1 7 .  (1  978). 

Jiusto,  J.r..,  "Considerations  in  the  Optical  Cnaractenzation  of  the 
Atmosphere,"  Atmos.  Sci.  Lab.,  White  Sands  Missile  Range,  NM, 

ASL-CR- 79-0 1  00-3  ,  (1979). 

Junge,  C.E.,  "Vertical  Profiles  of  Condensation  Nuclei  in  the  Stratosphere, 
Handbook  of  Geophysics  and  Space  Environments,  J.  Meteor.,  j_s,  p.  DU, 
(1961). 

Kampe,  H.J.A.,  and  Weikmann,  H  .  K .  ,  "Physics  of  Clouds,"  Meteor. 
Monographs,  3./  IP/  (1957). 

Kaveney,  W.J.  ,  Feddes ,  R.G.,  and  Liou ,  K-N,  "Statistical  Inference 
of  Cloud  Thickness  from  NCAA  4  Scanning  Radiometer  Data,"  Monthly- 
Weather  Review,  1 05 .  pp.  99-107  ,  (1977). 

Kessler,  R.  ,  "Kinimatical  relations  between  Wind  and  Precipitation 
Distributions  II,"  J.  Meteor.,  j_8,  (1961). 

Kessler,  E.  ,  "On  the  Distribution  and  continuity  of  water  substance  in 
Atmospheric  Circulations,"  Meteor.  Monogr.  ,  1 0 ,  84pp.,  (1969). 

Kessler,  E.  ,  "Model  of  Precipitation  and  Vertical  Air  Current,"  Tellus, 
XXVI,  pp.  519-542,  (1974). 

Kennedy,  J.S.,  and  Norberg,  W.  ,  "Circulation  Features  of  the  Stratospneie 
Derived  from  Radiometric  Temperature  Measurements  with  the  TIROS  VII 
Satellite,"  J.  Atmos.  Sci.,  24  ,  6,  (1967). 


98 


Kessler,  E.,  Ill,  and  Atlas,  D. , 
Aero/Space  Eng.,  Dec.,  (1959). 


Model  Precipitation  Distributions," 


Khrgian,  A.K.  and  Mazin,  l.P. 
Raindrop  Distribution  Spectra," 
lation)  17 .  36t46,  (1956). 


,  "Analysis  of  Methods  of  Characterizing 
Tr.  Tsents.  Aeroloa.  Observ.  Moscow,  (trans- 


King,  M.D.,  and  Herman,  B.M.,  "Determination  of  the  Ground  Albedo 
and  the  Index  of  Absorption  of  Atmospheric  Particulates  by  Remote 
Sensing.  Part  II:  Theory,"  J.  Atmos.  Sci.,36.,  (1979). 

Kneizys,  F.X.,  Shettle,  E.P.,  Gallery,  W.O.,  Chetwynd,  J.H.,  Jr., 

Abreu,  L.W.,  Selby,  J.E.A.,  Fenn,  R.W.,  and  McClatchley,  R.A.  , 
"Atmospheric  Transmittance/Radiance:  Computer  Code  LOWTRAN  5,"  Optical 
Physics  Division,  Hanscom  AFB,  AFGL-TR-80-0067 ,  (19B0). 

Knollenberg,  R.G.,  "Techniques  for  Probing  Cloud  Structures,"  IFAORS 
Abstract  Digest,  Clouds:  Their  Formation,  Optical  Properties  and  Effects, 
(1980). 

Kocmond,  W.C.,  Pilie,  R.J.,  Eadie,  W.J.  and  Mack.E.J.,  "Recent  Fog 
Modification  Research  at  Cornell  Aeronautical  Laboratory,"  Second  National 
Conference  on  Weather  Modification,  Santa  Barbara,  (1970). 

Kreitzberg,  C.W.,  and  Rasmussen,  R.G.,  "Treatment  of  Cloudiness  in 
Regional  Scale  Numerical  Weather  Prediction  Models,"  Preprint  Volume, 
Conf.  on  Cloud  Physics  and  Atmos.  Electricity,  Issaquah,  WA,  Amer. 
Meteor.  Soc.,  506-513,  (1978). 

Kreitzberg,  C.W.,  Mount,  W.D.,  and  Fow,  B.R.,  "Preliminary  Evaluation 
of  Meteorological  Models  for  Moisture  Depiction  and  Prediction  for  Electro- 
Optical  Applications,"  Final  Report  for  Atmos.  Sci.  Lab.,  White  Sands 
Missile  Range,  NM,  Contract  DAAG29-76-D-01 00  ,  (1979). 

Krishr.amurti ,  T.N.,  Ramanathan,  H-L  P.  ,  Pasch  R.J.,  Mohnari,  J., 
"Cumulus  Parameterization  and  Rainfall  Rates  I,"  Monthly  Weather  Review, 
108.  (1980). 


99 


Kumai,  M.  ,  "Arctic  Fog  Droplet  Size  Distribution  and  Its  Effect  on  Ligh‘ 
Attenuation,"  J.  Atmos.  Sci.,  3j3,  (1973). 

Kunkel,  B.A. ,  "Fog  Drop-Size  Distributions  Measured  with  a  Laser 
Hologram  Camera,"  J.  Appl.  Meteor.,  H),  (1 Q71  ) . 

Lala,  G.G.,  Mandel,  E.  ,  and  Jlusto,  J.E.,  "A  Numerical  Evaluation  of 
Radiation  Fog,"  J.  Atmos.  Sci.,  .32  ,  720-728,  (1975). 

Langille,  R.C.,  and  Gunn,  K.L.S.,  "Quantitative  Analysis  of  Vertical 
Structure  in  Precipitation,"  J.  Meteor.,  5.,  (1948). 

LaRocca,  A.].,  and  Turner,  R.E.,  “Atmospheric  Transmittance  and  itadianct 
Methods  of  Calculation,"  Contract  Nos.  N00014-73-A-0321-0002  and 
N  000 1 4-7  4-C- 0285  ,  Research  Ins.  of  Ml,  Ann  Arbor,  (1975). 

Laws,  J.O.,  and  Parsons,  D.A.  ,  "The  Relation  of  Raindrop-size  to 
Intensity,"  Transactions,  American  Geophysical  Union,  24.,  Pt.  2,  (1943). 

Lee,  I.Y.,  Hanel,  G.  ,  and  Pruppacher,  H.R.,  "A  Numerical  Determination 
of  the  Evolution  of  Cloud  Drop  Spectra  du*>  to  Condensation  on  Natural 
Aerosol  Particles,"  J  Atmos.  Sci.,  3_7,  (1980). 

Lewis,  R.E.,  " Meceorological  aspects  of  Aircraft  Icing,"  Compendium  of 
Meteor.  T.  F.  Malone,  ed.  AMS,  Boston,  (1951). 

Lowe,  P.R.,  "An  Approximating  Polynomial  for  Computation  of  Saturation 
Vapor  Pressure,"  J.  Appl.  Meteor.,  .1_6.,  p.  100-103  ,  (1977). 

Low,  D.H.R.,  "A  Theoretical  Investigation  of  Cloud/Fog  Extinction 
Coefficients  and  Their  Spectral  Correlations,"  Atmos.  Phys.,  52.,  1, 

(1979). 

Low,  R.D.H.,  "Fog  Evolution  in  the  Visible  and  Infrared  Spectral  Regions 
and  Its  Meaning  in  Optical  Modeling,"  Atmospheric  Sciences  Laboratory, 
White  Sands  Missle  Range,  NM,  (1979). 

Lovejoy,  S.,  and  Austin,  G.L.,  "The  Estimation  of  Rain  from  Satellite- 
Borne  Microwave  Radiometers,"  Q.J.R.M.S.,  1 06 .  pp.  255-276,  (1980). 


MacDonald,  J.E.,  "The  Saturation  Adjustment  in  Numerical  Modeling  of 
Fog,"  J.  Atmos.  Sci.,  20,  pp.  476-478,  (1963). 

Malm,  W.C.,  Walther,  E.G.,  and  Cudney,  R.A.,  "The  Effects  of  Water 
Vapor,  Ozone  and  Aerosol  on  Atmospheric  Turbidity,"  J.  Appl.  Meteor., 

16,  (1977). 

Matejka,  T.J.,  Houze,  R.A.  ,  Jr.,  and  Hobbs,  P.V.,  "Microphysics  and 
Dynamics  of  Clouds  associated  with  Mesoscale  Rainbands  in  Extratropical 
Cyclones,"  Q.J.R.M.S.,  106 .  pp.  29-56,  (1980). 

Mathai,  C.V.,  Harrison,  A.W.  ,  and  Mathews,  T.  ,  "Aerosol  Particle  Size 
Distribution  (0.1  -  1.0  m)  During  the  Chinooks  of  19  79  over  Calgary, 
Canada,"  J.  Appl.  Meteor.,  19,  (1980). 

Marshall,  J.S.,  and  Palmer,  W.M.,  "The  Distribution  of  Raindrops  with 
Size,"  J.  Meteor.,  _5,  pp.  165-166,  (1948). 

Mason,  B.J.,  D.Sc.,  F.R.S.,  "The  Physics  Of  Clouds Clarendon  Press, 
Oxford,  Second  Edition,  (1971).,  671  pp. 

McCartney,  E.J.,  "Optics  of  the  Atmosphere.  "John  Wiley  &  Sons,  New 
York,  407  pp. 

McClatchey,  R.A.,  and  Selby,  J.E.A.,  "Atmospheric  Attenuation  of  Laser 
Radiation  From  0.76  to  31.25  m,"  Optical  Physics  Lab.  Hanscom  AFB, 

AFCRL-TR-74-0003 ,  (1974). 

McClatchley,  R.A.,  Fenn,  R.W.,  Selby,  J.E.A.,  Volz,  F.E.,and  Garing, 
J.S.,  "Optical  Properties  of  the  Atmosphere"  (Third  Edition),  AF  CRL-72- 
0497,  1972. 

Medhurst,  R.G.,  "Rainfall  Attenuation  of  Centimeter  Waves:  Comparison  of 
Thoery  and  Measurement,"  Tefe  Transactions  on  antennas  and  Propagation, 
ed.  Falcone,  pp.  550-556,  (1955). 

Meyer,  M.B.,  Jiusto,  J.E.,  and  Lala,  G.G.,  "Measurements  of  Visual 
Range  and  Radiation-Fog  (Haze)  Microphysics,"  J.  Atmos.  Sci.,  37.,  PP* 
622-629,  (1980). 

Monohan,  H.H.,  "A  Summary  of  Meteorological  Data  Sources  in  the  United 
States  and  Studies  Relating  to  the  Climatology  of  Central  Europe,"  ASL 
Internal  Report,  (1979). 


101 


1 


Monohan,  H.  H.,  and  Cionco,  R.M.,  "An  interperative  Review  of  Exis'ir.g 
Capabilities  for  Measuring  and  Forecasting  Selected  Weather  Variables 
(Emphasizing  Remote  Means),  ASL  Report,  ASL-TR-0001  ,  (197«). 

Mooradian,  G.C.,  and  Geller,  M.,  "Pulsed  Blue-Green  Propagation  Through 
Clouds  and  Fogs",  IFAORS  Abstract  Digest,  Clouds:  Their  Formation, 

Optical  Properties  and  Effects,  (1980). 

Mount,  W.D.,  and  Fow ,  B.R.,  "Millimeter  and  microwave  obscuration  for 
10  mm  to  10  cm  Wavelengths.  Final  Rept.  GAC-70-003  to  Eatteile  Columbus 
Laboratories.  (1978). 

Mount,  W.D.,  Fow,  B.R.,  and  Mount,  B.D.,  "Cloud/Icing/Visibility 
Analysis  System  (CIVAS)  and  Forecast  Development  Plan,"  Final  Rept.  GAC- 
79-03-05-02  to  White  Sands  Missile  Range.  (1980). 

Mount,  W.D.,  Fow,  B.R.,  Mount,  B.D.,  and  Rapp,  J.O.  ,  "Sensitivity 
Study  of  CFAS  and  CFAR  Objective  Analysis  Techniques,"  Final  Rept.  GAC- 
78-CR-002  to  White  Sands  Missile  Range.  (1978). 

Muelle  ,  E.A.  and  Jones,  D.M.A.,  "Dropsize  distributions  in  Florida," 
Preprints  8th  Weather  Radar  Conf.  ,  San  Francisco,  Amer.  Meteor.  Soc .  , 

2‘»9-  .05  .  (1960). 

Nor  mi,  C.C.,  Mosher,  F.R.,  Hinton,  B.  ,  Martin,  D.W.,  Santek,  D.  ,  ant. 
Kuhlow,  W.  ,  "A  Model  for  Calculating  Desert  Aerosol  Turbidity  over  the 
Oceans  from  Geostat  ionary  Satellite  Data,"  J.  Appl.  Meteor.,  J_9 ,  6,  (l'c"0 

Ohtake,  T.  ,  "Observations  of  Size  Distributions  of  Hydrometers  Through  the 
Melting  Layer,"  J.  Atmos.  Sci.  ,  .26.,  (1969). 

Pasceri,  R.E.,  and  Friedlander,  S.K.,  "Measurements  of  the  Particle  Size 
Distribution  of  the  Atmospheric  Aerosol:  II.  Experimental  Results  and  Dis¬ 
cussion,  "  J.  Atmos.  Sci.,  22_,  (1965). 

Patterson,  E.M.,  Kiang,  C.S.,  Delany,  A.C.,  Wartburg,  A.F.,  Leslie, 
A.C.D.,  and  Huebert,  B.J.,  "Global  Measurement  of  Aeiosols  in  Remote 
Continental  and  Marine  Regions:  Concentrations,  Size  Distributions,  and 
Optica]  Properties,"  (paper  80C1078),  J.  Goephys.  Res.  .85,  C1  2,  (1  980). 


102 


Pasqualucci,  F.  ,  "Drop-Size  Distributions  in  Unicell,  Multicell,  and 
Squall  Line  Storms,"  Nat,  Conf.  on  Cloud  Physics,  6.2,8,  Boulder,  Colo., 
(1976). 

Passarelli,  R.E.JJr.,  "Notes  and  Correspondence:  Comments  on  "Size 
Distributions  of  Precipitation  Particles  in  Frontal  Clouds,"  J.  Atmos.  Sci., 

21,  (1980). 

Patterson,  E.M.,  and  Gillette,  D.A.,  "Measurements  of  Visibility  VS 
Mass-Concentration  for  Airborne  Soil  Particles,"  Atmos.  Environment, 

1 1  ,  pp.  193-196,  Pergamon  Press,  (1977). 

Patterson,  E.M.,  and  Gillette,  D.A.,  "Commonalities  in  Measured  Size 
Distributions  for  Aerosols  Having  a  Soil-Derived  Component,"  J.  Geophys. 

Sci .  ,  82,  1  5,  (1977). 

Pilie,  R.J.,  Mack,  E.J.,  Kocmond,  W.C.,  Eadie,  W.J.,  and  Rogers,  C.W., 
"The  Life  Cycle  of  Valley  Fog.  Part  II:  Fog  Microphysics,"  J.  Appl.  Meteor., 
Li,  (1  975). 

Pilie,  R.J.,  Mack,  E.J.,  Kocmond,  W.C.,  Rogers,  C.W.,  and  Eadie,  Vi.J., 
"The  Life  Cycle  of  Valley  Fog.  Part  1:  Micrometeorological  Characteristics," 

J.  Appl ,  Meteor.,  _M,  (1975). 

Pinmck,  R.G.,  Hoihjelle,  D.L.,  Fernandez,  G.,  Stenmark,  E.B.,  Lindberc, 
J.D.,  Hoidale,  G.B.,  and  Jennings,  S.G.,  "Vertical  Structure  in  Atmos¬ 
pheric  Fog  and  Haze  and  Its  Effects  on  Visible  and  Infrared  Extinction," 

T.  Atmos.  Sci.,  3_5  ,  (1978). 

Platt,  C.M.R.,  "Transmission  of  Submillimeter  Waves  Through  Water  Clouds 
and  Togs,"  J.  Atmos.  Sci.,  2_7,  (1970). 

Pratte,  F.J.,  Huffaker,  R.M.,  Lawrence,  T.R.,  and  Loveland,  R.  ,  "System 
Considerations  of  a  Long  Range  Coherent  Lidar  Wind  Sensor  With  Particular 
Emphasis  on  Real  Weather  Effects,"  NOAA  Tech.  Memo,  ERL  WPL-42,  Wave 
Propagation  Laboratory,  Boulder  Colorado,  (1979). 

Prospero,  J.M.,  and  Carlson,  T.N.,  "Vertical  and  Aerosol  Distribution  of 
Saharan  Dust  Over  the  Western  Equa*orial  North  Atlantic  Ocean,"  J.  Geo¬ 
phys.  Res.,  Jl,  p.  5255-5265,  (1972). 


103 


Prui-pacher ,  H.R.,  "The  Microstructure  of  Atmospheric  Clouds  an: 
tation,"  IFAORS  Abstract  Digest,  Clouds:  Their  Formation,  Ordeal 
and  Effects,  ( 1 9  30) . 


-Tecip.- 
Pr  o  per  *  1  e  s 


Pruppacher ,  H.R.,  and  Xlett ,  J.D.,  "Microphysics  of  Clouds  and  si  il¬ 
lation,  "  p.  9-27,  (1980), 


fj 


Pueschel,  R.F.,  and  Kuhn,  P.M.,  "Infrared  Absorption  of  Tropospheric 
Aerosols:  Urban  and  Rural  Aerosols  of  Phoenix,  Arizona,"  Coop:  vs. 
Res  .  ,  J30,  2  1  ,  (1975) . 


Cuenzel,  H .  ,  "Determination  of  Size  Distribution  of  atmospheric  Aerosol 
Particles  from  Spectral  Solar  Radiation  Measurements,  "  J.  Geophys.  Mcs .  , 
75 .  1  5,  (1970). 

Radke,  L.F.,  Hobbs,  P.V.,  and  Eltgroth,  M.W.,  "Scavenging  of  Acr  sm 
Particles  by  Precipitation,"  J.  Appl.  Meteor.  J_9,  (198Q). 

Reitan,  C.H.,  "Droplet  Size  Measurements  in  Convective  Clouds,  In 
■Artificial  Stimulation  of  Rain,  (eds .  H.  Weickmann  and  W,  Smith),  p. 
rergamon  Press,  NY,  (1957). 


Roach,  W.T.,  Brown,  R.  ,  Caugney,  S.J.,  Garland,  J.A.,  ar.d  Readings, 

C.J.,  "The  Physics  of  Radiation  Fog:  1  -  A  Field  Study.,  "  C.I..R.  M.S.  ,  '  .  ^ , 
432,  (1976). 

Roach,  W.T.,  "On  the  Effect  of  Radiative  Exchange  on  the  Growth  by 
Condensation  of  a  Cloud  or  Fog  Droplet,"  Q.J.R.M.S.,  1 02 .  432,  (1976). 

Rogers,  R.R.,  "A  Short  Course  in  Cloud  Physics."  2nd  Edition,  2 7  5  pp.  , 

Pergamon  Press,  NY,  (1979) 

Ryan,  B.F.,  "A  Model  to  Derive  Precipitation  Patterns  and  Vertical  Velo¬ 
cities  in  Large-Scale  Rain  Systems,"  Q.J.R.M.S.,  106,  pp.  201-212,  (19891. 


Savage,  R.C.,  "The  Radiative  Properties  of  Hydrometers  at  Microwave 
Frequencies,"  J.  Appl.  Meteor.,  1 7 .  p.  904-907,  (1978). 


Selby,  J.E.A.,  Shettle,  E.  P.  ,  and  McClatchey,  R.A.  ,  "Atmospheric  Trans¬ 
mittance  From  0.25  to  28.5  m:  Supplement  LOWTRAN  3B  (1976),"  AFGL-TR- 
76-0258  Environmental  Research  Papers,  Hanscom  AFB,  (1976). 


104 


Shettle,  E.P.,  and  Fenn,  R.W.,  "Models  for  the  Aerosols  of  the  Lower 
Atmosphere  and  the  Effects  of  Humidity  Variations  on  Their  Optical  Prop¬ 
erties,"  AFGL-TR-79-0214  Env .  Res.  Papers  No.  676,  Air  Force  Geopnys- 
ics  Laboratory.  (19  79). 

Snettle,  E.P.,  "Comment  on  Atmospheric  Aerosol  Size  Spectra:  Rapid 
Concentration  Fluctuations  and  Bimodaiity"  by  T.E.  Graedel  and  J.P. 
Praney,"  J.  Geophys.  Res.,  J30,  21  ,  (1975). 

Silverman,  B.A.,  and  Sprague,  E.D.,  "Airborne  Measurement  of  In-Cloud 
Vis'bility,”  p.  27  1-276  ,  Second  National  Conference  on  Weather  Mod¬ 
ification,  Santa  Barbara,  (1970). 

Silverman,  B.A.,  "An  Eulerian  Model  of  Warm  Fog  Modification,"  Air 
Force  Cambridge  Res.  Labs.,  Bedford,  MA. 

Smith,  P.L.,  Myers,  C.G.,  and  Orville,  H.D.,  "Radar  Reflectivity  Factor 
Calculations  in  Numerical  Cloud  Models  Using  Bulk  Parameterization  of 
Precipitation,").  Appl.  Meteor.  .14,'  (1975). 

Spilhaus ,  A.F.,  "Drop  Size,  Intensity,  and  Radar  Echo  of  Rain,"  J.  Me’., 
3,  lcl,  (1948). 

Srivastava,  R.C.,  "Size  Distribution  of  Raindrops  Generated  By  Their 
Breakup  and  Coalescence,"  J.  Atmos.  Sci.  ,  28 .  p.  401-415,  (1971). 

Srivastava,  R.C.,  "Parameterization  of  Raindrop  Size  Distributions,"  J. 
Atmos.  Sci.,  35  f  108-1  17  ,  (1978). 

Srivastava,  R.C.,  and  Atlas,  D.  ,  "Growth,  Motion  and  Concentration  of 
Precipitation  Particles  in  Convective  Storms,"  J.  Atmos.  Sci.,  2_6,  (1969). 

Steinvall,  O.  ,  "Computed  Me  Scattering  Properties  for  Laser  Wavelengths 
in  Various  Atmospheric  Media,"  FCA-2-C-2662-B1 -B3 ,  Res.  Inst,  of  Nat. 
Def.  41  p.  ,  (1971) . 

Tag,  P.M.,  Johnson,  D.B.,  ENS,  Hindman,  E.E.  II,  "Engineering  Fog- 
Modification  Experiments  by  Computer  Modelling,"  Second  National  Conf. 
on  Weather  Modification,  Santa  Barbara,  (1970). 


105 


Toon,  O.B.,  and  Pollock,  J.B.,  "A  Global  Average  Model  of  Atm  m  nr: ic 
Aerosols  tor  Radiative  Transfer  Calculations,”  J,  Appl.  Me’eor.,  i  5 . 
p.  225-246,  (1976). 

Trusty,  G.L.,  and  Cosdcn,  T.H.,  “An  Aerosol  Measuremen*  System  For 
Laser/Aerosol  Interaction  Studies,"  Proceedings  of  the  Op?icai-Submilhrr<*tef 
Atmospheric  Propagation  Conference,  Vol.  1,  (1976). 

Twomey ,  S.,  and  Warner,  J.  ,  "Comparison  of  Measurements  of  Cloud 
Droplets  and  Cloud  Nuclei,”  J.  Atmos.  Sci.  ,  p.  402-403,  6, 

Twomey,  S.,  "Aerosol  Size  Distribution  by  Multiple  Filter  Measurements  ,  ” 

J.  Atmos.  Sci.,  p.  1073-1079,  3_3_ ,  (1976). 

Valley,  S.L.,  "Handbook  of  Geophysics  and  Space  Environments,  "  Air 
Force  Cambridge  Research  Laboratories,  Chapt.  5,  p.  5-1  to  5-27,  (1965). 

Varley,  D.I.,  Lt.  Col.,  USAF,  "Microphysical  Properties  of  a  Large 
Scale  Cloud  System,  1-3  March  1978,"  AFGL-TR-80-0002 ,  Environmental 
Research  Papers  No.  690,  Air  Force  Geophysics  Laboratory,  (1980). 

Wang,  P.K.,  and  Iruppacher,  H.R.,  "Acceleration  to  Terminal  Velocity  'f 
Cloud  and  Raindrops,"  J.  Appl.  Me'eor.,  p.  275-280,  JL_6 ,  (1  977). 

Warner,  J.  ,  and  Twomey,  S.,  "The  Production  of  Cloud  Nuclei  by  Cane 
Fires  and  the  Effect  on  Cloud  Droplet  Concentration,"  J.  Atmos.  Sci.,  24 , 

6 ,  (19  6  7). 

Warner,  J.  ,  "The  Microstructure  of  Cumulus  Cloud.  Part  I.  General 
Features  G  the  Droplet  Spectrum,"  J.  Atmos.  Sci.,  p.  1049-1059,  26 . 
(1969). 

Warner,  J.,  "Notes  and  Correspondence:  On  the  Spatial  Variability  of  Drop¬ 
lets  in  Cloud,"  J.  Atmos.  Sci.,  d.  541-545  ,  3_6<  (1979). 

Weibel,  G.E.,  and  Dressel,  H.O.,  "Propagation  Studies  in  Millimeter- Wave 
Link  Systems,"  Proceedings  of  the  IEEE,  55,  4,  (1967). 

Wexler ,  R.  and  Atlas,  D.  ,  "Moisture  Supply  and  Growth  of  Stratiform  Pre¬ 
cipitation,"  J.  Meteor.,  15,  p.  531-539  ,  (1958). 


106 


Yamamoto,  G.  et  al.,  'Table  of  Scattering  Function  of  Infrared  Radiation 
for  Water  Clouds,"  NOAA  Technical  Report  NESS  57  Cat.  No.  COM-71- 
503  12.  GPO,  Washington,  DC,  (1971). 

Young,  K.C.,  "The  Evolution  Of  Drop  Spectra  Due  To  Condensation, 
Coalescence,  and  Breakup,"  J.  Atmos.  Sci.,  965-973,  32,  (1975). 

Zissis,  G.J.,  and  Wolfe,  W.L.,  et  al.,  "The  Infrared  Handbook,"  En¬ 
vironmental  Research  Institute  of  Michigan,  (1978). 

Squires.P.  : 

-The  Microstructure  and  Colloidal  Stability  of  Warm  Clouds.  Pt.  I.  The 
Relation  Between  Structure  and  Stability.,"  Tellus  JJ3,  (195B). 

-The  Microstructure  and  Colloidal  Stability  of  Warm  Clouds.  Pt.ll  The 
Causes  of  the  Variations  in  Microstructure,"  Tellus,  _1_0,  (1958). 

-The  Microstructure  of  Cumuli  in  Maritime  and  Continental  Air,"  Tellus,  8, 
443,  (1956). 

-The  Spatial  Variation  of  Liquid  Water  and  Droplet  Concentration  In  Cumuli," 
Tellus,  H),  372,  (1959). 

L'.S.  Standard  Atmos pnere,  1976,  p.  44-47  ,  NOAA,  NASA,  USAF,  Washington, 
DC,  1976.  Superintendent  of  Documents. 


"Penetrability  of  Haze,  Fog,  Clouds  and  Precipitation  by  Radiant  Energy 
over  the  Spectral  Range  0.1  Micron  to  10  Centimeters,"  CNA  Naval  Warfare 
Analysis  Group,  Study  61,  (1969). 

Aerosol  Measurement,  1979  ,  Univ.  Presses  of  Florida.  D.A.Lundgren  et  al. 

Inadvertent  Climate  Modification.  Report  on  the  Study  of  Man’s  Impact  on 
Climate  (SM1C)  by  M.I.T.,  1971,  M.I.T.  Press. 

Gunn,  R.  and  Kinzer,  D.  D.  ,  "The  Terminal  Velocity  of  Fall  of  Water  Droplets 
In  Stagnant  Air,"  J.  Afeteorol.  6,  p.  243,  (1949). 

aufm  Kampe,  H.J.,  and  Weickmann,  H.K.,  "Physics  of  Clouds,"  Meteorological 
Monographs,  3.,  p.  182,  (1957). 

Weikmann ,  H.K.,  "Physics  of  Precipitation,"  Meteorological  Monographs,  3, 

P.  266,  (1957). 


107 


